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SECTION  I 


INTRODUCTION 

No  critical  review  exists  of  the  thermal  decomposition 
mechanisms  of  polynitroaromatic  explosives.  This  void  is 
surprising  in  view  of  the  fact  that  2,4,6-trinitrotoluene  (TNT) 
has  been  available  as  a  pure  material  since  about  1870  and  is  the 
most  widely  used  explosive.  For  any  explosive  the  thermal 
decomposition  mechanisms  and  resultant  products  are  important 
details  because  they  affect  the  opportunity  to  predict 
systematically  the  feasibility  of  large-scale  synthesis,  the  long¬ 
term  stability  for  purposes  of  storage,  and  the  sensitivity  to 
various  stimuli,  such  as  heat  and  mechanical  impact.  Although 
safety  and  performance  of  explosives  are  intimately  tied  to  the 
thermal  decomposition  process,  knowledge  of  the  mechanisms  does 
not  automatically  enable  improvements  in  safety  and  performance  to 
be  made.  However,  without  mechanistic  information,  the  behavior 
of  explosives  remains  highly  empirical. 

Robertson  and  Yoffe  (Reference  1)  and  Wenograd  (Reference  2) 
were  among  the  earliest  to  connect  thermal  decomposition  processes 
to  initiation  phenomena.  Moreover,  the  thermal  decomposition 
mechanisms  and  products  are  believed  to  influence  the  performance 
of  the  explosive,  such  as  the  Chapman- Jouguet  (C-J)  pressure  (the 
pressure  at  the  rear  of  the  reaction  zone  of  the  detonation  front) 
and  the  detonation  velocity.  For  example  Cook,  al .  (References 
3  and  4)  proposed  that  the  initial  reaction  kinetics  were 
important  in  determining  the  detonation  velocity.  Zeman  and  co¬ 
workers  expanded  this  notion  by  attempting  to  correlate  the 
thermal  decomposition  kinetics  of  polynitroaromatic  explosives  at 
lower  temperature  with  molecular  structure  (Reference  5) ,  heat  of 
explosion  (Reference  6) ,  detonation  characteristics  (References  5 
and  7),  and  thermal  stability  (Reference  8).  The  critical 
temperature  of  a  thermal  explosion  is  related  to  Arrhenius 
kinetics  by  the  Frank-Kamenetskii  model  (References  9-12) . 

A  very  large  number  of  polynitrobenzene  derivatives  are 
known,  of  course.  All  reasonable  trinitrobenzene  derivatives  will 
explode  if  subjected  to  an  energy  input  exceeding  a  threshold 
value.  TNT  achieves  its  special  position  partly  because  of  its 
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compatibility  with  other  materials  and  low  hygroscopic ity,  and 
partly  because  -of  its  good  thermal  stability,  relatively  low 
sensitivity  to  impact  and  friction,  and  high  explosion  power. 

These  characteristics  influenced  the  German  Army  to  adopt  TNT  in 
place  of  picric  acid  (PA) ,  which  is  more  hydroscopic,  chemically 
reactive,  and  sensitive  to  impact.  Other  nitrobenzene  derivatives 
having  practical  value  as  explosives  are  triamino-trinitrobenzene 
(TATB) ,  2,4, 6-tr initrophenylmethy lnitramine  (Tetry 1 ) , 
hexanitrostilbene  (HNS),  and  tetranitrodibenzo-l,3a,4,6a- 
tetrazapentalene  (TACOT) ,  which  is  used  in  two  isomeric  forms,  Z 
and  T,  that  differ  in  the  position  of  the  -N02  groups  on  the 
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the  rings.  TATB  has  excellent  thermal  stability  and  is  especially 
noteworthy  for  its  extreme  insensitivity  to  impact.  Likewise, 
TACOT  has  exceptional  thermal  stability.  It  is  able  to  retain 
predictable  explosive  characteristics  after  ten  minutes  at  350*C 
making  it  useful  for  deep  well  oil  exploration.  Meyer  (Reference 
13)  is  a  good  source  of  practical  aspects  of  explosives. 

A  large  amount  of  work  has  been  conducted  on  decomposition 
processes  of  polynitroaromatics .  Of  the  many  conclusions  reached, 
some  are  contradictory,  many  are  perplexing,  some  are  in 
unavailable  government  documents,  and  others  are  in  untranslated, 
poorly  circulated  (at  least  in  the  West)  Russian  literature.  This 
article  is  an  attempt  to  assess  critically  the  literature  of 
decomposition  processes  of  polynitroaromatic  explosives  and  to 
come  to  some  conclusions  about  the  dominant  thermal  mechanisms  in 
the  bulk  phase. 

First,  it  is  necessary  to  outline  why  thermal  decomposition 
mechanisms  of  explosives  have  produced  divergent  literature. 
Treacheries  of  the  materials  aside,  the  complications  and 
realities  of  thermolysis  processes  of  relatively  large  molecules 
in  their  bulk  state  and  the  restrictions  on  measurements  are 
highlighted  next  in  Sections  II.  1-4. 
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SECTION  II 


SPECIAL  CONSIDERATIONS  OF  EXPLOSIVE  DECOMPOSITION  CHEMISTRY 

The  results  and,  ultimately,  the  description  of  thermal 
decomposition  of  an  explosive  depend  on  the  amount  and 
configuration  of  the  sample,  the  rate  and  magnitude  of  energy 
loading,  the  diagnostics  employed,  and  the  interpretative  scheme. 
Experimental  conditions  are  frequently  dramatically  different  from 
one  laboratory  to  another  which  complicates  comparison  of  the 
data. 

1.  Decomposition  Depends  on  the  Conditions  of  Thermal  Stimulus 

The  numerous  decomposition  reactions  of  a  complex  organic 
molecule  in  the  bulk  condensed  phase  each  have  their  own  rate 
constant  and  sensitivity  to  temperature.  The  temperature  to  which 
the  bulk  material  is  heated  and  the  heating  rate  controls  the  time 
that  the  products  have  to  form.  At  lower  temperatures,  processes 
with  lower  activation  energies  and  greater  exothermicity  are 
favored.  Time  is  available  for  stable  products  to  grow  in 
concentration  via  the  complex  reaction  network  and  to  serve  as 
catalysts.  At  higher  temperatures,  the  higher  activation  energy 
and  more  endothermic  processes  increase  in  importance.  However, 
the  products  have  a  shorter  time  to  form  and  may  •  :  may  not  be  the 
same  as  products  formed  at  lower  temperature.  At  least  the 
product  distribution  will  be  different.  Consequently, 
descriptions  of  decomposition  processes  depend  significantly  on 
the  temperature  and  the  heating  rate,  a  fact  that  has  led  to  many 
contradictory  assertions.  Aside  from  the  chemistry,  the  physical 
details  about  the  material,  such  as  its  thermal  conductivity,  heat 
capacity,  hardness,  void  volume,  void  size,  etc.,  contribute  to 
the  details  of  the  decomposition  process,  especially  during 
mechanical  stress.  The  conversion  of  mechanical  energy  into 
chemical  effects  (Reference  14)  is  being  addressed  at  an 
increasingly  sophisticated  level  (References  15-19) . 

a.  Low  Heating  Rate/Low  Temperature 

Programmed  heating  over  minutes  or  hours  to  a  chosen 
temperature  while  monitoring  the  change  of  heat  flow  or  pressure 
is  the  basis  of  differential  scanning  calorimetry  (DSC) , 
differential  thermal  analysis  (DTA) ,  and  manometry.  The 
decomposition  rates  extracted  tend  to  apply  to  the  lower 
activation  energy  processes  and  reflect  the  overall  (global) 
process,  rather  than  a  particular  step.  Zeman  and  co-workers 
(References  5-8)  claim  that  correlations  exist  between  the 
activation  energies  during  the  induction  period  before 
autocatalysis  takes  over  and  detonation  properties  pertaining  to 
high  pressure  and  temperature,  such  as  the  heat  of  explosion 
(Reference  6) ,  C-J  pressure,  and  detonation  velocity  (Reference 
7) .  They  suggest  that  the  rate  of  early  decomposition  process  is 
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extendable  to  high  temperature  end  pressure  conditions.  The 
products  that  are  isolated  by  solvent  extraction  of  TNT  heated  to 
a  temperature  characteristic  of  slow  decomposition  are  numerous 
and  frequently  difficult  to  identify  (References  20  and  21) . 

b.  Moderate  Heating  Rate/High  Temperatures 

Heating  times  in  the  range  of  20  -  500  Msec  to  a  chosen  high 
temperature  can  be  achieved  by  T-jump  methods  (References  2,22). 
Wenograd  (Reference  2)  quotes  heating  rates  of  an  explosive  packed 
into  a  hypodermic  needle  at  2.5  x  10f  K/sbc.  These  rates  are 
undoubtedly  optimistically  fast  because  models  indicate  that  the 
rate  of  heat  transfer  to  the  explosive  cannot  keep  pace  with  the 
heat  that  is  believed  to  be  supplied  (Reference  23} .  Heatinr 
times  of  <1000  Msec  to  a  chosen  temperature  are  important  be  a 
they  simulate  conditions  that  exist  during  mechanical  impact  t  a 
during  combustion.  In  fact,  a  frequently  used  index  of  the 
sensitivity  of  an  explosive  to  mechanical  impact  is  the  drop 
weight  impact  test.  The  impact  sensitivity  is  the  height,  H30, 
that  a  weight  of  given  mass  must  be  dropped  onto  the  sample  to 
produce  an  explosion  50%  of  the  time.  The  duration  of  heating 
during  impact  is  200  -  300  Msec  (Reference  24)  and  the  pressure 
achieved  is  7  -  15  kbar.  In  addition  to  the  chemical  reactions, 
hso  is  affected  by  a  score  of  other  variables  including  the 
particle  size,  hardness  of  the  crystals,  the  defect  number,  the 
defect  volume,  the  crystal  size  distribution,  the  atmospheric 
conditions,  and  the  operator  technique.  Given  the  many  factors 
that  influence  the  impact  sensitivity,  numerous  shots  are  averaged 
while  holding  constant  as  many  of  the  variables  as  possible.  Not 
suzprisingly,  some  scatter  in  the  data  occurs. 

Robertson  and  Yoffe  (Reference  1)  concluded  that  chemical 
decomposition  during  impact  is  very  similar  to  that  at  less  rapid 
heating  rates.  This  notion  was  extended  by  Wenograd  (Reference  2) 
who  found  that  the  impact  sensitivity  of  a  series  of  diverse 
explosives  roughly  correlated  with  the  temperature  to  which  the 
sample  could  be  raised  in  250  Msec.  Thereby,  he  too  concluded 
that  impact  sensitivity  was  linked  to  thermally  induced  chemistry. 
The  temperature  achieved  by  several  explosives,  including  TNT,  in 
this  time  exceeded  1000 *C.  It  should  be  emphasized  that  this 
result  does  not  link  thermal  stability  with  impact  sensitivity, 
because  thermal  stability  is  controlled  by  lower  activation  energy 
reaction  branches  which  become  de-emphasized  in  favor  of  higher 
activation  energy  branches  at  higher  temperatures. 

The  drop  weight  impact  test  does  not  provide  chemical  details 
about  the  thermal  processes.  To  achieve  this  information,  a 
technique  such  as  T-jump/FTIR  spectroscopy  (Reference  22)  can  be 
used.  This  technique  is  able  to  separate  and  detect  some  of  the 
dominant  reaction  branches  during  rapid  heating,  and  to  determine 
the  temperature  dependence  of  some  of  these  branches  (References 
22,25) . 
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c.  High  Heating  Rate/High  Tempera ture/High  Pressure 

The  use  of  a  shock  wave  to  heat  a  material  affords  a  heating 
time  to  a  high  temperature  of  less  than  1  nsec .  For  example,  the 
chemical  peak  corresponding  to  95%  decomposition  TNT  in  a 
stationary  detonation  has  a  duration  of  0.4  nsec  (Reference  26). 
Pressures  for  shocked  condensed  phase  materials  are  in  the  range 
of  30  -  20r>  kbar.  The  shock  sensitivity  of  many,  but  not  all, 
polynitro&'omatics  for  which  data  are  available  follows  the  same 
trend  as  the  impact  sensitivity  (Reference  27)  suggesting  that  the 
shock  initiated  chemistry  is  similar  to  impact  initiated 
chemistry.  However,  a  considerable  body  of  Russian  literature 
indicates  that  this  is  not  always  the  case  (Reference  28) . 

An  up-to-date  overview  of  shock  and  detonation  physics  of 
explosives  is  available  elsewhere  (Reference  29) .  In  practice 
shock  compression  of  a  solid  or  liquid  can  be  induced  in  a  number 
of  ways  that  usually  employ  a  flyer  plate.  The  most  common  method 
is  to  use  a  standard  explosive  to  transmit  a  shock  wave  through  a 
buffer  plate  to  the  explosive  of  interest.  Tests  known  as  NwedgeN 
and  "gap"  are  well-known  to  workers  in  the  field.  For  example, 
the  Navy  small  scale  gap  test  uses  a  booster  explosive  in  a 
standard  configuration  to  transmit  a  high  pressure  shock  wave 
through  a  Plexiglas  attenuator  to  the  explosive  packed  at  a 
particular  density  compared  to  the  theoretical  maximum  density. 

The  shock  sensitivity  is  determined  by  the  pressure  at  which  an 
explosion  occurs  50%  of  the  time.  For  example,  P90  would  be  the 
initiation  threshold  pressure  for  an  explosive  packed  at  90% 
theoretical  maximum  density.  A  comparison  of  various  detonation 
sensitivity  tests  is  given  by  Urizar,  fit  al.  (Reference  30) . 
Electric  spark  sensitivity  has  also  been  compared  to  shock  and 
impact  sensitivity  (Reference  31) . 

For  studies  of  unimolecular  reactions  in  the  gas  phase  at  low 
pressure,  a  laser  pulse  (Reference  32)  or  a  shock  wave  created  by 
a  burst  of  gas  pressure  (References  33,34)  can  be  used  to  heat 
molecules  to  a  high  temperature  in  <1  Msec.  Temperatures  in 
excess  of  800 *C  can  be  achieved  by  shock  tube  and  laser  pyrolysis. 
Lower  and  higher  activation  energy  processes  can  be  expected  to 
compete  effectively  with  one  another  under  these  conditions. 

In  summary,  the  reaction  temperature  and  the  heating  rate 
play  an  important  role  in  how  the  thermal  decomposition  process  of 
a  complex  molecule  occurs.  Some  reactions  can  be  expected  to  be 
present  over  a  wide  range  of  temperatures  while  others  will  not. 
Consequently,  considerable  care  is  needed  when  comparing  the 
merits  of  reaction  pathways  derived  from  measurements  made  under 
diverse  conditions. 
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2.  Experimental  Diagnostics  for  Chemistry  of  Explosives 

Interpretations  of  the  chemistry  that  is  important  in  an 
explosive  depend  on  the  diagnostic  method,  because  the  time 
resolution,  the  pressure,  the  temperature,  the  sample  quantity, 
and  the  phase  are  major  variables.  The  relatively  high  vapor 
pressure  of  many  nitroaromatic  compounds  makes  them  amenable  to 
study  by  gas  phase  as  well  as  condensed  phase  analytical 
techniques . 

a.  Dynamic  Analysis 

Ideally  the  course  of  the  reactions  could  be  followed 
directly  in  real  time.  When  this  condition  is  achieved  for  an 
explosive  (Reference  35) ,  the  spectral  vinJow  is  narrow  and  the 
amount  of  detailed  information  is  limited.  More  commonly, 
simulation  of  conditions  is  required  to  optimize  the  chemical 
detail. 

Manometry  is  used  extensively  in  Russia,  especially  by 
Maksimov,  to  determine  the  rate  of  gas  evolution  as  a  function  of 
temperature  for  many  nitroaromatic  compounds.  Typically,  a  15  - 
20  cm3  glass  vessel  fitted  with  a  Bourdon-type  compensation 
manometer  is  immersed  in  a  thermostated  bath.  Evaporation  of  the 
compound  is  sought  so  that  the  decomposition  rate  can  be  measured 
in  the  gas  phase.  However,  transfer  of  heat  from  the  glass 
surface  to  the  nitroaromatic  molecule  at  the  temperatures  used  is 
primarily  by  conduction.  The  molecule  must  collide  with  the  wall 
where  heterogeneous  catalysis  can  be  appreciable.  Nevertheless,  a 
large  percentage  of  the  rate  data  for  aromatic  nitrocompounds  has 
been  obtained  by  this  method. 

DSC  and  DTA  have  been  applied  extensively  to  obtain  rate  data 
for  the  condensed  phase  decomposition  of  energetic  nitroaromatic 
compounds.  These  techniques  extract  the  rate  of  the  processes 
from  the  overall  heat  flow  at  relatively  low  temperatures  and 
heating  rates.  Details  of  the  kinetics  are  sometimes  over¬ 
shadowed  by  the  fact  that  energetic  materials  can  generate  heat 
rapidly  and  the  instrument  is  limited  by  its  ability  to  detect 
heat  changes  instantly.  As  a  result,  the  rate  data  are  sensitive 
to  the  sample  size.  As  with  manometry,  the  global  reaction  rate 
is  measured.  The  reactions  that  are  responsible  for  the  rate  must 
be  inferred. 

Mass  spectrometry  is  more  broadly  applicable  to  the  study  of 
thermal  decomposition  mechanisms  of  nitroaromatic  compounds  than 
it  is  for  many  other  energetic  materials.  The  foundation  of  this 
notion  lies  in  the  extensive  work  of  Fields  and  Meyerson 
(References  36,37),  who  investigated  the  mass  spectra  of  numerous 
nitroaromatic  compounds  and  concluded  that  the  early  fragmentation 
patterns  closely  resemble  the  thermal  decomposition  routes.  This 
finding  can  not  be  applied  broadly  to  other  energetic  materials. 
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With  nitramines,  for  example,  the  species  produced  by  thermal 
routes  and  those  generated  by  electron  impact  must  be  carefully 
differentiated  (Reference  38) . 

Electron  impact  mass  spectrometry  yields  only  the  fragments 
and  not  their  mechanism  of  formation  (Reference  39) .  Tandon  mass 
spectroscopy  (NS/MS)  used  in  conjunction  with  collision  induced 
dissociation  (CID)  elucidates  the  fragmentation  pathways  of 
individual  ions.  The  principle  is  that  mass  selected  primary  ions 
are  fragmented  in  a  reaction  cell  as  a  result  of  collisions  with 
an  inert  gas.  These  daughter  ions  are  then  analyzed  by  a  second 
MS,  thereby  enabling  determination  of  the  individual  fragmentation 
pathways  of  the  parent  ions  (Reference  40) . 

Electron  spin  resonance  potentially  yields  mechanistic 
information,  but  it  is  sensitive  only  to  radicals  that  form  and 
degrade  during  thermal  decomposition  in  the  condensed  phase.  Very 
low  concentrations  (<  10'*  M)  are  readily  detected.  There  is  no 
way  to  prove  that  the  radicals  observed  are  the  key  species  in  the 
decomposition  mechanism.  Non-radical  and  ionic  pathways  are 
always  possible.  Assignment  of  structure  and  composition  of  the 
species  is  sporting,  especially  for  complicated  condensed  phase 
decomposition  reactions  that  produce  closely  related  radicals. 

b.  Post  Reaction  Analysis 

Rapid-scan  FTIR  spectroscopy  is  a  near  real-time  analytical 
technique  with  a  wide  spectral  window.  The  gaseous  decomposition 
products  at  sub-atmospheric  and  super-atmospheric  (-100  atm)  can 
be  determined  at  heating  rates  up  to  2  x  10**C/sec  (Reference  22) . 
Structure/decomposition  relationships  can  be  developed  (Reference 
41)  and  the  sequence  of  product  formation  outlines  the 
decomposition  mechanism  (References  22,25).  The  products  detected 
are  not  the  elementary  radicals,  but  the  more  thermally  and 
kinetically  stable  species. 

Unimolecular  decomposition  in  the  gas  phase  can  be  induced  at 
high  temperature  by  use  of  a  pulsed  laser  (References  32,42)  or  a 
shock  impulse  (References  33,34).  The  quenched  products  can  be 
analyzed  by  a  variety  of  techniques,  such  as  by  GC-MS. 

Essentially  the  full  range  of  analytical  techniques  can  be 
applied  to  analyzing  species  in  a  quenched  sample  that  have  been 
generated  by  heating  or  mechanically  abusing  the  material.  The 
caveat  is  that  the  most  reactive  intermediates  which  are  central 
to  the  decomposition  mechanism  at  elevated  temperature  will 
rearrange  and  react  further  upon  cooling  and  during  time  delayed 
analysis.  The  thermodynamically  and  kinetically  more  stable 
species  therefore  are  always  detected.  These  species  are  only  a 
historical  record  of  the  most  relevant  dynamic  processes. 


7 


3.  Complexity  of  Thermal  Decomposition  in  the  Condensed  Phase 

Unperturbed  unimolecular  decomposition  kinetics  can  be 
determined  only  for  molecules  in  the  gas  phase  at  low  pressure. 
Materials  that  are  useful  as  explosives  or  propellants  almost 
always  begin  to  react  in  the  condensed  phase.  Thus,  while  gas 
phase  studies  of  individual  molecules  can  be  useful  guides  to  the 
behavior  in  other  phases,  the  most  important  reactions  arise  when 
the  molecules  are  in  the  condensed  phase  in  the  presence  of  like 
molecules.  For  instance,  by  using  the  oxidation  of  hydrocarbons 
as  an  analogy,  the  liquid  phase  reaction  is  rarely  unimolecular, 
but  the  propagation  steps  can  be  similar  in  the  liquid  and  the  gas 
phase  (Reference  43) . 

Qualitative  accounts  of  the  decomposition  process  of  bulk  TNT 
and  TATB  indicate  that  determination  of  the  mechanisms  of 
polynitroaromatic  derivatives  in  the  presence  of  neighboring 
molecules  is  an  extraordinary  difficult  problem.  Robertson 
(Reference  44)  described  the  behavior  of  TNT  when  heated  to  a 
constant  temperature  in  a  closed  ampule.  TNT  melted  at  81*C  and 
partially  volatilized,  followed  by  a  quiescent  period  with  no  gas 
evolution.  Next,  gas  evolution  began  accompanied  by  darkening 
followed  by  more  rapid  gas  evolution  and  then  explosion.  During 
the  period  of  darkening,  the  amount  of  gas  is  small  (about  1  mole 
per  mole  of  TNT)  and  catalysts  form  (Reference  20) .  Catalano  and 
Crawford  (Reference  45)  observed  that  exothermic  decomposition  of 
TATB  begins  at  about  302 *C  and  found  that  the  reaction  stage  in 
the  342  -  357 *C  range  necessitated  the  use  of  five  different  rate 
expressions.  These  observations  are  stark  reminders  that  the 
aggregation  of  many  time-dependent  chemical  and  physical  processes 
occurs  during  the  thermal  decomposition  of  a  material  in  the 
condensed  phase.  Despite  the  complications,  major  temperature 
segments  of  the  decomposition  regime  follow  Arrhenius  behavior, 
equation  1.  E,  is  the  activation  energy,  A  is  the  frequency 
factor,  and  k  is  the  reaction  rate. 

k(T)  *  Ae'SJFT  <x> 

The  major  perturbations  to  the  elementary  decomposition 
reactions  of  an  isolated  polynitrobenzene  molecule  when  the 
molecule  is  part  of  a  bulk  condensed  phase  are: 

(1)  Bimolecular  reactions  (References  46,47)  and 
intermolecular  associations  (Reference  48) ; 

(2)  Homogeneous  and  heterogenous  (References  32,49,50) 
catalysis  (References  20,34,44,51)  by  reaction  products; 

(3)  Heterogeneous  catalysis  by  the  supporting  surface  and  the 
reactor  walls  (References  7,32,44,48,49,52-56). 
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Intermolecular  association  and  atom  transfer  have  been 
proposed  by  Maksimov  (Reference  48)  to  be  responsible  for  the  fact 
that  the  rate  of  decomposition  of  TNT  in  the  melt  phase  is  about 
ten  times  faster  than  in  the  gas  phase  or  in  solution.  The 
reverse  is  true  of  most  other  energetic  molecules,  such  as  nitrate 
esters,  hexahydro-1, 3, 5-trinitro-s-triazine  (RDX) ,  and  octahydro- 
1,3,5,7-tetranitrotetrazacine  (HMX)  (Reference  48).  In  fact,  the 
specific  rate  constants  for  many  condensed  phase  decomposition 
reactions  of  polynitroaromatics  are  10s  -  io10  sec'1,  which  is  below 
the  limit  of  10u  sec'1  for  nany  unimolecular  reactions  and  is 
suggestive  of  bimolecular  processes  (Reference  48) . 

An  important  factor  contributing  to  the  faster  rate  of 
condensed  phase  reactions  compared  to  the  gas  phase  is  catalysis 
by  the  reaction  products.  For  this  reason  virtually  every  kinetic 
study  of  decomposition  of  polynitroaromatic  explosives  in  the 
condensed  phase  since  an  early  report  (Reference  57)  has  found 
that  an  autocatalytic  expression  is  required  to  fit  the  rate  data. 
Of  course,  very  small  amounts  of  catalyst  can  markedly  influence 
the  rate  of  most  reactions.  In  fact,  Maksimov  (Reference  58) 
proposed  that  impurities  in  the  original  sample  sharply  change  the 
rate  of  decomposition  of  nitroaromatic  compounds.  Irrespective  of 
the  impurities,  many  compounds  (References  20,44,51,59,60) 
including  essentially  any  source  of  H*  or  other  radicals, 
accelerate  the  rate  of  decomposition  of  TNT  (Reference  51) .  TNT 
and  most  of  its  decomposition  products  produce  radicals  at  the 
10*4  M  level  (Reference  61)  when  heated  below  300*C.  Several 
specific  decomposition  products  of  TNT  have  been  tested  and  are 
found  to  catalyze  decomposition  of  TNT.  These  include  2,4,6- 
trinitrobenzaldehyde  (Reference  20),  2 , 4-dinitroanthranil  (2,4- 
DNAn)  (References  51,59),  an  ill-defined  explosive  "coke"  from 
long-term  thermal  degradation  (Reference  20) ,  and  elemental  carbon 
(Reference  62) .  The  role  of  catalysis  by  low  levels  of 
decomposition  products  is  not  readily  separable  from  the  role  of 
intermolecular  associations  (Reference  48)  mentioned  above. 
Therefore,  catalysis  rather  than  intermolecular  associations  may 
be  responsible  for  the  faster  rate  of  decomposition  of  TNT  in  the 
condensed  phase  compared  to  the  gas  phase. 

The  influence  of  catalysis  by  the  walls  of  the  reactor  and 
the  supporting  surface  of  the  sample  has  been  mentioned  by  many 
authors  as  a  potential  factor  in  the  reaction  rate  (References 
7,32,44,48,49,52-54,63),  but  has  been  investigated  systematically 
in  only  a  few  studies.  For  example,  Maksimov  (Reference  53)  found 
that  increasing  the  surface  area  of  the  glass  reaction  vessel  with 
the  aid  of  capillaries  greatly  accelerated  the  initial  rate  of 
decomposition  of  trinitrobenzene  derivatives  in  the  vapor  phase. 
Since  the  wall  and  bulk  phase  decomposition  mechanisms  are  thought 
to  be  the  same  (References  46,53,63),  it  is  difficult  to  separate 
the  role  of  catalysis  by  the  decomposition  products  from  the  role 
of  surface  catalysis. 
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On*  of  the  more  helpful  empirical  insights  for  sorting  out 
thermal  decomposition  schemes  of  trinitrobenzene  derivatives  in 
the  condensed  phase  is  the  observation  that  initial  fragmentation 
patterns  in  the  El  mass  spectrum  seem  to  parallel  thermal 
decomposition  schemes  in  the  gas  phase  (References  36,37,64,65). 
Moreover,  it  is  helpful  that  the  early  thermal  decomposition 
products  found  in  the  gas  phase  closely  resemble  the  stable 
intermediates  formed  by  heating  the  condensed  phase.  While  the 
elementary  steps  leading  to  these  molecules  are  undoubtedly 
altered  somewhat  by  the  phase,  the  similarity  of  the  stable 
species  in  the  gas  phase  and  the  stable  intermediates  in  the 
condensed  phase  supports  a  role  for  them  in  the  decomposition 
sequence,  at  least  under  certain  conditions.  However,  because  of 
catalysis  by  the  decomposition  products  and  intermolecular 
interactions  as  noted  above,  there  is  no  reason  to  assume  that  the 
reaction  details  of  a  nitrobenzene  derivative  in  the  gas  phase  or 
in  dilute  solution  will  be  the  same  as  the  bulk  material. 

4.  Uncertainty  about  Determining  Molecular  Control  over 

Explosive  Behavior 

It  was  noted  in  Section  I  that  there  are  strong  reasons  to 
believe  that  thermochemistry  is  in  some  way  connected  to  impact 
and  shock  initiation  of  explosives  (References  1-12,27,66).  The 
desire  to  understand  the  controlling  molecular  mechanisms  of 
impact  and  shock  sensitivity  of  an  explosive  has  led  to  many 
published  correlations  of  explosive  data  with  molecular 
properties.  The  most  extensive  of  these  compilations  for 
nitroaromatic  compounds  are  by  Kamlet  and  Adolph  (Reference  24), 
Zeman  (References  5-8,54),  and  Storm,  £&  al.  (Reference  27).  A 
number  of  other  correlations  of  molecular  and  crystal  parameters 
with  explosives  behavior  have  been  noted  for  the  limited,  but 
attractively  related,  set  of  the  four  compounds  TNB,  MATB,  DATB 
and  TATB.  However,  this  series  is  also  an  excellent  illustration 


MATB  DATB  TATB 

j 

of  the  problems  of  correlating  of  molecular  variables  with 
explosive  properties  in  an  attempt  to  find  the  controlling 
mechanism.  For  TNB,  MATB,  DATB  and  TATB  these  variables  are 
categorized  according  to  whether  the  basis  is  predominantly  the 
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explosive  characteristics  (Table  1) ,  the  composition  and  crystal 
properties  (Table  2)  or  the  molecular  and  electronic  structure 
(Table  3) .  Most  of  the  entries  are  self-explanatory ,  but  a  few 
require  brief  description.  The  impact  sensitivity,  H50,  the  shock 
initiation  threshold  pressure,  P90,  and  the  Chapman- Jouguet 
pressure,  P&,,  were  defined  in  Sections  I  and  II.  1.  AH,,  is  the 
heat  of  explosion.  In  Table  2,  pN  is  the  density  of  -N02  groups 
in  the  crystal  lattice  defined  by  equation  2  (Reference  72)  where 

PAT*-—  (2) 


Z  is  the  number  of  molecules  in  the  unit  cell,  M  is  the  molecular 


100  (2N0-NH-2Ne-2Neon) 


weight,  and  V  is  the  unit  cell  volume.  The  oxidant  balance,  OB100, 
is  defined  according  to  equation  3  as  the  number  of  equivalents  of 
oxygen  in  100  gms  of  a  CHNO  explosive  required  to  convert  all  of 
the  H  to  H20  and  C  to  CO.  N  is  the  number  of  atoms  of  each 
element  in  the  formula.  The  melting  point  is  partly  a  reflection 
of  the  intermolecular  cohesive  forces.  In  Table  3,  the 
electrostatic  potential  calculated  at  the  mid-point  of  the  C-NOz 
bond,  V>id,  determined  by  a  Mulliken  orbital  population  analysis  at 
the  level  of  CINDO  (Reference  73)  and  ab  initio  SCF-MO  with  STO-3G 
(Reference  74)  or  STO-5G  (Reference  75)  basis  sets.  The  number  of 
aromatic  ir  electrons  was  calculated  at  the  STO-3G  level  (Reference 
67).  The  shake-up  promotion  energy  separation,  ASPE,  of  N(ls)  of 
the  -N02  groups  was  obtained  by  x-ray  photoelectron  spectroscopy 
(Reference  76) .  The  average  bond  distances  were  obtained  from  the 
crystal  structures  (Reference  68-71) . 

Perusal  of  the  columns  of  Tables  1-3  and  of  the  Tables 
relative  to  one  another  reveals  that  a  trend  exists  in  all 
columns.  To  emphasize  the  similarity,  the  correlation 
coefficient,  R,  (Reference  2)  and  the  confidence  level  of 
statistical  significance  according  to  the  t  test  (Reference  77) 
equation  (4),  are  given  in  Table  4  for  all  columns  plotted  versus 
the  impact  sensitivity,  H50.  To  exceed  95%  confidence,  the  value 
of  t  for  four  data  points  (n)  must  exceed  3.18.  In  all  cases 


except  one,  this  confidence  level  is  achieved.  In  fact,  twelve 
data  sets  exceed  99%  confidence  (t  >  5.84).  The  only  case  that 
fails  95%  exceeds  90%  confidence  (t  >  2.35).  Furthermore,  two 
sets  of  data  that  are  linearly  related  to  a  third  are  linearly 
related  to  each  other.  Therefore,  among  the  18  sets  of  data  in 
Tables  1-3,  there  are  153  positive  correlations.  A  few  of  course 


11 


Table  1.  Sensitivity  and  Performance  Parameters 


Compound 

H,0,  cm 

P,o,  Jcbar* 

Pa,  GPa<* 

AH*,.,  kJ  g 

TNB 

100* 

15 

22.9 

4.734 

MATB 

177* 

28 

24.7 

4.263 

DATB 

320* 

46 

26.7 

4.138 

TATB 

490* 

70 

29.3 

3.990 

a.  Reference  24 

b.  Estimated,  but  frequently  used,  see  References  27  and  72. 

c.  Reference  5 

d.  Reference  7 

e.  Reference  6 


Compound 

Table  2.  Compositional  and  Crystal  Parameters 

MW,  g  mol*1  p,  g  cm'3  pN*,  g  cm"3  AH£,cal  g*1  mp,  K 

OBioo 

TNB 

213 

1.688* 

0.653 

l 

a 

H 

• 

OS 

396 

-1.46 

MATB 

228 

1.760b 

0.646 

-78 

461 

-1.75 

DATB 

243 

1.837® 

0.628 

-122.6 

559 

-2.06 

TATB 

254 

1.938d 

0.623 

-143 

6231 

-2.33 

a.  Reference  68 

b.  Reference  69 

c.  Reference  70 

d.  Reference  71 

e.  Reference  72 

f.  Not  well  characterized. 
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Molecular  and  Electronic  Paraneters 


Table  4.  The  Coefficients  of  Correlation  and  the  Confidence  Level 
based  on  the  Nt  test"  for  Plots  of  all  Entries  in  Tables 


1-3  versus 

the  Impact  Sensitivity  Index,  Hso. 

Parameter 

Correl.  Coeff,  R2 

t  (equation 

P»o 

0.998 

31.6 

Pcj 

0.992 

15.75 

AH.X 

0.773 

2.61 

W 

0.952 

6.30 

P 

0.991 

14.83 

pN 

0.927 

5.04 

AH*f 

0.931 

5.20 

mp 

0.972 

8.33 

OB10o 

0.971 

8.18 

STO-5G 

0.982 

10.45 

V^,,,  3TO-3G 

0.972 

8.49 

V^,  CINDO 

0.919 

4.76 

ASPE 

0.950 

6.17 

No.  w  e* 

0.968 

7.78 

^C-»02 

0.981 

10.17 

^H-0 

0.955 

4.55 

<?c-e 

0.994 

18.2 

a.  Significant  at  confidence  level  90%  >  2.35?  95%  >  3.18;  99%  > 
5.84. 
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are  dependent  or,  at  least,  interconnected.  The  density  and  C-J 
pressure  are  related  through  the  detonation  velocity,  D,  by 
equation  5.  The  molecular  weight  enters  into  OB100  by  equation  3. 

Pcj  •  0 . 25pZ>3  (5) 

The  increase  in  the  average  C-C  and  N-0  distances  is  paralleled  by 
the  decrease  in  the  average  C-N02  distance  as  -NH2  groups  are 
added.  In  some  cases  the  correlations  are  pure  folly  even  though 
R  (Reference  2)  is  high.  For  example,  melting  point  vs.  shake-up 
promotion  energy  separation,  heat  of  formation  vs.  density, 
molecular  weight  vs.  the  C-N02  mid-point  electrostatic  charge  are 
among  the  more  chemically  ridiculous  correlations. 

On  the  other  hand,  several  of  the  correlations  of  the  columns 
in  Tables  1-3  have  been  proposed  to  identify  crystal  and  molecular 
factors  that  influence  the  impact  and  shock  sensitivity  of  TNB, 
MATB,  DATB  and  TATB.  The  correlation  of  Hso  in  Table  1  with  the 
shake-up  promotion  energy  separation  (References  76,78)  (Table  3) 
was  proposed  as  evidence  of  the  importance  of  excited  states  in 
determining  the  impact  sensitivity  (Reference  78) .  Correlation  of 
Hso  with  V^d  in  Table  3  has  been  suggested  to  show  that  the  C-N02 
bond  strength  and  homolysis  play  a  role  (References  74,75)  or  are 
rate-controlling  (Reference  73)  in  the  impact  and  shock 
sensitivity  of  the  TNB,  MATB,  DATB,  TATB  series.  Table  4  reveals 
that  the  C-N02  bond  distance  correlates  as  well  or  better  with  H50 
than  V.!,,.  In  fact,  all  of  the  trends  in  Table  3  are  principally 
the  result  of  systematic  ground  state  electronic  and  structural 
trends  along  this  series.  That  is,  the  C-C,  C-NOz  and  N-o  bond 
distances  are  no  more  or  less  "controlling"  of  the  sensitivity 
than  any  other  factor  at  this  level  of  analysis.  Systematic 
differences  in  the  bond  distances  figure  importantly  into  the 
computed  electrostatic  potentials  and  measured  relative  orbital 
energies . 

Other  correlations  of  data  in  Table  2  with  H50  have  been 
made.  Odiot  (Reference  72)  noted  that  the  density  of  -N02  groups 
in  the  crystal  correlated  with  Hs0.  Kamlet  and  Adolph  (Reference 
24)  correlated  OB100  with  the  log  H50  for  TNB,  MATB,  DATB  and  TATB 
and  other  compounds.  The  trends  implied  that  the  oxidizing 
potential  of  the  molecule  may  be  important.  Table  3  reveals  that 
for  this  series  of  compounds,  Hso  correlates  as  well  with  OB100  as 
does  log  H50.  Storm,  sfe  al.  (Reference  27)  noted  that  within 
Table  3,  Hso  and  P90  correlate  well.  Rogers,  ££  al.  (Reference  67) 
proposed  that  the  increasing  amount  of  intermolecular  hydrogen 
bonding  as  -NH2  groups  are  added  is  responsible  for  the  decrease 
in  shock  sensitivity.  As  intermolecular  hydrogen  bonds  increase 
in  number,  they  absorb  energy  from  the  shock  front  and  prevent  it 
from  being  localized  in  the  molecule  (Reference  67) .  Although 
difficult  to  quantify  with  the  data  available,  a  qualitative 
indication  of  the  trend  in  intermolecular  association  in  the 
series  is  the  melting  point,  and,  possibly,  the  density.  Both 
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increase  as  -NH2  groups  are  added  to  the  ring  and  the  shock 
sensitivity,  P,0,  decreases. 

An  inescapable  conclusion  of  Tables  1-4  is  that  single  column 
correlations  of  the  molecular,  lattice,  and  explosive  properties 
for  TNB,  MATB,  DATB  and  TATB  are  not  especially  helpful  for 
identifying  the  primary  chemical  mechanisms  of  shock  and  impact 
sensitivity.  Correlations  within  a  given  Table  are  the  soundest. 
It  is  especially  risky  to  pick  a  column  of  one  table  and  correlate 
it  with  a  column  in  another  table.  For  example,  any  column  of 
Table  3  is  as  justifiable  as  any  other  to  correlate  with  Hso  in 
Table  1.  As  a  result,  previous  conclusions  about  molecular 
features  that  control  shock  and  impact  sensitivity  in  this  series 
of  compounds  are  contradictory:  excited  electronic  states 
(Reference  78)  vs.  C-N02  electrostatic  potential  (References  73- 
75)  vs.  intermolecular  hydrogen  bonding  (Reference  67)  vs.  oxidant 
balance  (Reference  24) .  With  153  positive  correlations  among  the 
data,  single  correlations  are  unlikely  to  provide  the  true 
factors.  Perhaps  all  of  these  factors  (or  none  of  them) 
contribute  to  the  sensitivity  to  some  degree. 

Given  the  enormous  barriers  to  obtaining  meaningful  data  on 
the  condensed  phase  decomposition  mechanisms  of  nitroaromatic 
explosives  (temperature  effects,  pressure  effects,  shortcomings  of 
diagnostics,  catalysis,  bimolecular  processes,  rapid  heat 
generation,  parallel  properties,  et  .) ,  it  is  not  surprising  that 
contradicting  data  and  interpretations  are  common.  In  the 
following  sections  the  merits  of  kir.etics  and  mechanistic  models 
of  thermal  decomposition  are  evaluated. 


16 


SECTION  III 


GLOBAL  KINETICS  OF  THERMAL  DECOMPOSITION 

Most  knowledge  about  the  thermal  decomposition  kinetics  of 
explosives  is  obtained  from  rate  measurements  at  relatively  low 
temperature  (<  400 *C).  Kinetic  data  for  TNT  are  cited  in  this 
section  along  with  some  discussion  of  the  merits  of  measurement 
methods.  These  kinetic  data  are  used  in  later  sections  as  an  aid 
in  determining  the  reactions  mechanisms. 

As  mentioned  in  Section  II. 3 ,  the  decomposition  process  of 
bulk  TNT  involves  a  solid-to-liquid  phase  change.  Reactions  in 
the  liquid  phase  lead  to  a  small  amount  of  gas  and  non-volatile 
intermediate  products,  followed  by  other  reactions  that  produce  a 
large  quantity  of  lower  molecular  weight  gases.  It  is  widely 
agreed  (References  44,79,80)  that  the  overall  process  in  the 
liquid  phase  is  modeled  by  a  first-order  rate  expression  (6) , 

ln(l-a)  =  k(T)  t  (*) 

where  a  is  the  amount  of  sample  decomposed  and  t  is  time. 
Autocatalysis  and  self -heating  complicate  the  interpretation  of 
the  kinetics.  Therefore,  attempts  to  measure  k(T)  are  usually 
made  under  adiabatic  or  isothermal  conditions  whose  relationship 
to  k(T)  for  a  first-order  process  is  given  by  the  heat  balance 
equation  (7)  (Reference  9) .  X  Is  the  thermal  diffusivity,  c  is 
the  specific  heat,  and  Q  is  the  heat  of  reaction.  The  first  term 
of  equation  (7)  accounts  for  heat  loss  to  the  surroundings.  The 
second  term  is  the  heat  accumulated  by  the  explosive.  The  sum  of 
these  two  terms  equals  the  heat  generated  by  the  thermal 
decomposition.  Equation  (7)  can  not  be  solved  under  non- 

-XV2r  +  pC(dTfdt)  =  pQA(l-a)e~B‘,RT  <7> 

isothermal  non-adiabatic  conditions.  However,  under  adiabatic 
conditions  (X^^  *  O) ,  such  as  time-to-explosion  or  time-to- 
exotherm  (TTX)  experiments  (References  11,44,81)  the  apparent 
activation  energy  can  be  calculated  from  equation  (8) ,  where  t.  is 


mt.  •  -§&  n  <•) 

the  time-to-exotherm  at  temperature  T.  Under  isothermal 
conditions  (dT/dt  =  O) ,  such  as  is  achieved  by  evolved  gas  and  DSC 
methods,  Arrhenius  constants  are  obtained  from  equation  (1) . 

Figure  1,  which  is  an  idealized  representation  of  the  DSC 
thermogram  of  TNT,  can  be  subdivided  into  three  regions  (Reference 
82).  Based  on  equation  (7),  the  induction  period  represents  the 
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ACCELERATORY  DECAY  PHASE 
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Figure  1.  The  Energy  Release  Trace  for  Slow  Decomposition  of  a 
Nitroaromatic  at  a  Relatively  Low  Temperature. 


length  of  tine  that  the  sanple  is  thermally  quiescent  afterheating 
toa  set  temperature.  Activation  energy  determined  during  this 
period  may  have  the  best  chance  of  being  related  to  the  structure 
of  the  parent  molecule,  since  it  reflects  the  earliest  reaction 
events.  The  acceleratory  stage  is  the  time  during  which  catalysis 
and  self-heating  dominate  the  rate  of  decomposition.  The  decay 
phase  is  the  time  when  later  stage  decomposition  products  react. 

In  reality,  the  kinetics  of  decomposition  of  a  bulk  sample  of  an 
explosive  are  also  sensitive  to  the  mass  and  geometry  of  the 
sample.  One  reason  is  that  the  rate  of  heat  transfer  into  and  out 
of  the  sample  can  not  keep  pace  with  the  chemical  reaction. 

Hence,  the  heat  changes  detected  are  partly  a  measure  of  the 
reaction  rate  and  partly  a  measure  of  the  heat  transfer  rate. 

Mass  transport  may  also  occur  at  a  different  rate  than  the 
reaction  rate.  Under  these  circumstances,  agreement  among 
measured  rate  parameters  obtained  by  different  techniques  is 
unlikely.  Moreover,  association  of  the  rate  of  the  process  during 
any  of  the  phases  in  Figure  1  with  a  specific  reaction  is 
necessarily  tentative. 

The  Arrhenius  parameters  for  TNT  are  compiled  in  Table  5,  and 
fall  into  three  groups  according  to  the  stages  shown  in  Figure  1. 
These  data  apply  below  400 *C  and  are  not  necessarily 
representative  of  the  chemistry  at  high  temperature  which  is 
initiated  by  shock  and  impact.  Considerations  of  mechanisms  at 
higher  temperature  are  discussed  in  more  detail  in  Section  IV. 

Most  publications  in  the  past  15  years  mention  the  stage  of  Figure 
1  that  is  believed  to  be  studied.  If  not  noted,  the  data  were 
grouped  according  to  the  Ea  value.  Older  literature  rarely 
discusses  the  reaction  stage.  The  Arrhenius  constants  of  the 
induction  phase  tend  to  be  higher  than  the  other  stages  perhaps 
because  they  are  most  closely  related  to  a  bond  breaking  step. 
Maksimov  (Reference  84)  attributes  his  data  in  this  range  to 
unimolecular  processes.  In  the  acceleratory  stage,  the  Arrhenius 
constants  are  smaller  perhaps  because  of  catalysis  by  the 
decomposition  products.  Maksimov  (Reference  84)  attributes  his 
rate  data  in  this  range  to  bimolecular  processes.  The  measurement 
by  Urbanskii  and  Rychter  (Reference  85)  gives  the  lowest  E.  value 
of  all  probably  because  the  temperature  exceeds  the  boiling  point 
of  TNT  (345*C)  (Reference  44).  In  this  case,  heat  of  evaporation 
is  mixed  with  the  activation  energy.  The  value  of  Brill  and  James 
(Reference  87)  was  obtained  at  a  high  heating  rate  and  under  30 
atm  Ar  to  suppress  evaporation.  Wenograd's  low  value  of  Ea  was 
determined  at  high  temperature  and  high  heating  rate  on  a  semi- 
confined  sample.  It  is  nearly  impossible  to  achieve  sufficiently 
rapid  heat  transfer  and  temperature  sensing  under  Wenograd's 
conditions  to  obtain  an  accurate  value  of  Ea  (Reference  23).  In 
fact,  global  activation  energies  determined  at  high  heating  rates 
and  high  temperature  almost  always  tend  to  be  smaller  than  those 
at  lower  heating  rate  and  lower  temperature  (References  88,89). 

The  one  measurement  (Reference  82)  of  Arrhenius  parameters  that  is 
clearly  stated  to  be  of  the  decay  phase  is  similar  to,  but 
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Table  5.  Arrhenius  Pars— tsrs  for  Thar— X  Decomposition  of 
TUT.  Data  are  Arranged  Qualitatively  According 
to  the  Stage  of  the  Process  Measured. 

Induction  Phase 


T,  *C 

Ea,  keal/mol 

A,  sec'1 

Method* 

Ref 

203  -  283 

40.9  ±  1.6 

- 

ESR 

80 

298  *  490 

41.4 

1013.2 

TTX 

81 

237.5  -  276. 

5  43.4  ±  0.8 

^q12.19  t  0.32 

EG 

79 

140  -  180 

46  ±  3.7 

- 

EG 

83 

220  -  260b 

46.5 

1012* 

EG 

84 

245  -  269 

46.5  ±  1.5 

mm 

IDSC 

82 

240  -  255 

50.7  ±  2.4 

- 

DTA 

55 

Acccl.8Eat.ara ,  Phase 

390  -  450 

14 

- 

TTX 

85 

500  -  720 

18 

io‘-‘ 

TTX 

2 

-300 

27 

- 

TTX 

57 

203  -  283 

30.2  ±  0.6 

mm 

ESR 

80 

275  -  310 

32.0  ±  2.5 

- 

TTX 

44 

275  -  310 

34.4  1  2.5 

- 

EG 

44 

280  -  320 

34.5 

- 

EG 

53 

190  -  270 

34.6 

- 

EG 

86 

220  -  260b 

35.6 

10’ 

EG 

84 

390  -  407 

37.4 

- 

TTX 

87 

Decay  Phase 

245  -  269 

29.4  ±  1.4 

io’-’ 

IDSC 

82 

a.  EG  -  Manone trie  measurement  of  evolved  gas;  TTX  *  time-to- 
explosion  (exotherm) ;  IDSC  -  isothermal  differential 
scanning  calorimetry;  DTA  *  differential  thermal  analysis. 

b.  However,  rate  data  in  Figures  of  reference  84  are  given  at 
280-330‘C. 
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phase  is  similar  to,  but  slightly  less  than,  the  accelerator^ 
phase. 


The  methods  of  measurement  of  the  Arrhenius  data  in  Table  5 
reflect  all  of  the  complications  of  kinetics  determinations  of 
thermal  decomposition  of  bulk  phase  explosives.  Nanometric 
measurement  of  the  rate  of  gas  evolution  from  nitrobenzene 
explosives  has  been  used  extensively  in  Russia  especially  by 
Maksimov  and  co-workers.  They  described  manometry  as  measuring 
decomposition  in  the  vapor  phase,  but  noted  that  substituted 
nitrobenzene  compounds  occasionally  form  a  brown  residue  on  the 
reactor  wall  (Reference  53).  As  noted  in  Section  II. 3.,  catalysis 
by  this  residue  contributes  to  the  reaction  rate.  Use  of  the 
manometric  method  for  rate  measurements  has  been  critized 
(References  20,80-82)  because  the  amount  of  gas  generated  by 
trinitrobenzene  derivatives  depends  on  the  amount  of  non-volatile 
residue  that  forms  and  the  mechanism  of  the  reaction.  Moreover, 
many  substituted  trinitrobenzene  explosives  have  significant  vapor 
pressure  at  elevated  temperature  (Reference  90) .  The  vapor 
pressure  equation  for  TNT  is  given  by  equation  (9)  (References 
44,91).  Fortunately,  knowledge  of  the  vapor  pressure  enables  the 
sample  size  to  be  chosen  to  ensure  that  all  of  the  decomposition 

log  P{mm  Hg )  *  9.11-3850/T  (9) 

occurs  in  the  vapor  phase.  Table  5  reveals  that  the  activation 
energies  from  manometry  are  similar  to  those  obtained  by  other 
techniques.  Consequently,  criticisms  of  manometry,  while  probably 
valid,  do  not  disqualify  the  data  and  do  not  strengthen  confidence 
in  any  other  method  of  kinetics  measurement. 

Various  methods  of  measuring  the  time-to-exotherm  (explosion) 
involve  sensing  the  time  of  violent  release  of  heat  or  explosion. 
Equation  (8)  is  used  to  determine  the  apparent  activation  energy. 
Numerous  studies  employ  this  method  and  a  characteristic  of  most 
data  is  curvature  (References  2,80,81,92-95)  in  the  plot  of  lnt. 
vs.  1/T.  The  inability  to  control  the  heat  transfer  rate  and  to 
know  the  true  temperature  makes  the  value  of  Ea  sensitive  to  the 
assumed  temperature.  Perhaps  the  best-controlled  method  for 
measuring  the  length  of  the  induction  period  is  isothermal  DSC. 

Measurement  of  the  rate  of  growth  of  the  ESR  signal  of  a 
paramagnetic  species  has  been  used  to  determine  the  rate  of  the 
induction  and  acceleratory  phase  of  TNT  decomposition  (Reference 
80) .  Following  the  concentration  of  a  particular  species  in  this 
way  is  sounder  than  EG  or  DSC  measurements  of  the  global  rate,  but 
assumes  that  the  reaction  involving  the  paramagnetic  species 
occurs  at  the  same  rate  as  the  overall  decomposition  process. 

In  sum,  there  is  no  completely  unambiguous  method  to  measure 
the  rate  of  decomposition  of  an  energetic  material  in  the  bulk 
state.  Fortunately,  rough  agreement  exists  among  measurements  on 
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TUT  ltdt  toy  savaral  diffarant  aathoda.  Rata  data  axist  for 
dacoapoaition  of  othar  nitroaroaatic  compounds  (Rafarancaa  5,6,8, 
47,50,52, 53-56,84 , 86 , 96-106)  soaa  of  which  ara  usad  latar  in  this 
raport. 
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SECTION  IV 


EARLY  THERMAL  DECOMPOSITION  REACTIONS 

Several  of  the  thermal  decomposition  reactions  of 
polynitroaromatic  compounds  result  from  simple  bond  scission  or 
rearrangements  of  the  substituents  on  the  ring.  These  are  among 
the  first  few  reactions  and  thereby  initiate  all  subsequent  steps 
in  the  complex  decomposition  scheme.  This  section  evaluates  these 
reactions  beginning  with  the  nitro  group  and  concluding  with 
reactions  of  other  substituents  on  the  benzene  ring,  particularly 
those  ortho  to  the  -N02  group.  A  notable  feature  of  these 
reactions  is  the  fact  that  the  aromatic  ring  remains  intact. 
Destruction  of  the  aromatic  ring  is  a  later  stage  process 
described  in  Section  VI. 

1.  C-NOa  Homolysis 

Simple  homolysis  of  the  C-N02  bond  to  form  an  N02  radical  and 
an  aromatic  radical  by  equation  (10)  is  a  logical  reaction  on  the 


basis  of  the  importance  of  the  analogous  reaction  in 
nitroaliphatics,  nitramines  (R2NN02) ,  and  nitrate  esters  (RON02) 
(Reference  107) .  For  these  latter  compounds  the  X-N02  bond  is 
frequently  the  weakest  bond  in  the  molecule.  Two  observations 
suggest,  however,  that  simple  extension  of  this  analogy  to 
nitroaromatic  compounds  requires  qualification.  First,  the  bond 
dissociation  energy  (BDE)  of  nitrobenzene  derivatives  is  known  to 
be  about  71  kcal/mol  (Reference  32)  and  is  altered  by  only  1-4 
kcal/mol  by  adding  an  a-CH3  group  (Reference  32) ,  an  additional 
— N02  group  (References  32,52),  or  Cl,  Br,  I,  NHZ,  and  OH, 

(Reference  52)  to  the  ring.  The  C6H5-N02  BDE  is  significantly 
larger  than  that  of  an  aliphatic  carbon-bound  -N02  group  (about  60 
kcal/mol) ,  a  nitramine  (about  50  kcal/mol) ,  and  a  nitrate  ester 
(about  40  kcal/  iol)  (Reference  108).  Second,  although  N02(g)  has 
been  reported  from  decomposition  of  nitrobenzene  at  275 *c 
(Reference  109)  and  from  pyrolyzed  1, 3 -dinitrobenzene  (Reference 
110),  1,4-dinitrobenzene  (Reference  110),  and  1,3,5- 
trinitrobenzene  (Reference  111)  by  MS,  N02(g)  is  rarely  observed 
from  decomposition  of  substituted  nitrobenzene  compounds  in  the 
bulk  state.  This  behavior  sharply  contrasts  with  the  behavior  of 
nitroaliphatics,  nitramines,  and  nitrate  esters,  where  N02(g)  is  a 
readily  detected  product  (Reference  107) .  Both  of  these 
observations  suggest  that  nitroaromatic  compounds  may  have 
alternate  decomposition  channels  whose  activation  energies  lie 
below  that  of  C-N02  homolysis,  especially  when  other  substituents 
are  present  on  the  ring. 
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Table  6  lists  Arrhenius  parameter*  for  C-N02  homolysis  from 
unsubstitutsd  nitroaroaatic  compounds  (aquation  10)  in  the  gas 
phase.  The  process  is  first-order  (Reference  52).  Earlier 
literature  values  of  Ea  for  thermal  decomposition  of  nitrobenzene 
of  51-52  kcal/mol  (References  113,114)  appear  to  be  in  error.  The 
values  of  Ea  in  Table  6  are  very  similar  to  the  BDE  of  C6H5-H02, 
which  is  good  evidence  that  the  activation  energy  corresponds  to 
C-NO2  hoaolysis  (Reference  115) .  As  noted  above,  pyrolysis  of 
1,3-  and  1,4-dinitrobenzene  (Reference  110),  and  1,3,5- 
trinitrobenzene  (Reference  111)  in  the  gas  phase  produces  ions 
corresponding  to  the  successive  loss  of  N02.  The  rate  constants 
for  decomposition  of  these  three  compounds  in  the  gas  phase  at 
330*C  are  in  the  ratio  1:20:34  (Reference  53).  The  primary  loss 
of  N02  from  these  molecules  leaves  the  aromatic  fragment  with  the 
vibrational  and  electronic  internal  energy  to  drive  additional 
reactions  (Reference  111) . 

C— N02  homolysis  by  equation  (11)  is  a  dominant  decomposition 
channel  of  3-  and  4-nitrotoluene,  and  also  occurs  for  2- 
nitrotoluene,  (References  32-34,36,52,96,116).  Table  7  list  the 
Arrhenius  data  for  equation  (11).  While  only  rough  agreement 
exists  among  these  data,  the  Ea  values  are  similar  to  the  BDE  of 
C-N02.  The  relatively  large  A  factor  is  characteristic  of 


unimolecular  bond  scission  and  produces  a  positive  entropy  of 
activation.  Using  Maksimov's  data  (Reference  96)  As*  is  about  20 
cal/deg  for  3-NT  and  14.5  cal/deg  for  4-NT.  C-N02  homolysis  has 

been  suggested  to  be  initial  reaction  of  para-substituted 
nitrobenzene  compounds  on  the  basis  of  DSC  and  accelerating  rate 
calorimeter  data  where  the  temperature  at  the  onset  of 
decomposition  roughly  correlates  with  the  electron  withdrawing  and 
donating  effect  of  the  substitutent  (Reference  116a) . 

In  the  gas  phase  C-N02  homolysis  in  1,2 -NT  completes  with 
other  reactions.  These  reactions  have  been  addressed  by  Tsang  and 
co-workers  (References  33,34),  and  McMillen  and  Golden  and  co¬ 
workers  (References  32,42).  Equation  (12),  the  formation  of 
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Table  6.  Arrhenius  Data  for  C-NOa  Homolysis  of  Mono-,  Di-  and 
Trinitrobanzana  compounds 


Compound 

T,  • 

C 

Ea,  keal/mol 

A,  sac'1 

Raf. 

C^HsHOa 

410 

- 

480 

69.7 

±  1.4 

^q17.3  ±  0.3 

52 

455 

- 

530 

67.3 

±  4.7 

- 

112 

797 

- 

907 

65.7 

IqXS.3 

33 

827 

- 

977 

68.2 

±  1.7 

1015.5  t  0.5 

32 

l,3-(N02)aC2H* 

420 

- 

480 

68  1 

3 

10H.8  ±  0.5 

52 

827 

- 

977 

70  ± 

2.1 

101*.5  ±  0.5 

32 

1,4-  (NQa)  2C6H4 

420 

- 

470 

68.6 

±  2.5 

1017.1  ±  0.5 

52 

1,3,5-(N02)3C6H, 

380 

- 

470 

67.3 

±  0.8 

^q!7.2  ±  0.2 

52 
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Table  7.  Arrhenius  Data  for  C-NOa  Hoaolysis  of  Hononitrotoluene 
(NT)  and  Dinitrotoluene  (DNT)  Derivatives 


Coq^ound 

T, 

•c 

• 

Ea#  kcal/nol 

A,  sec'1 

Ref 

2 -NT 

797 

- 

907 

61.47 

10l*  81 

33 

827 

- 

977 

67.0  ±  2.2 

^q16.4  *  0.6 

32 

3 -NT 

300 

- 

350 

65 

1017  88 

96 

400 

- 

470 

68 

10“  ’ 

52 

4 -NT 

300 

- 

350 

61.5 

10“  73 

96 

400 

- 

470 

65.9 

10“7 

52 

797 

- 

907 

64.29 

10l*  M 

33 

827 

- 

977 

68.2  ±  2.0 

1014.9  ±  0.5 

32 

2, 4 -DNT 

827 

- 

977 

67.4  ±  1.7 

^q15.3  t  0.4 

32 
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anthranil,  was  found  to  be  an  important  reaction  channel  in  laser 
assisted  homogeneous  pyrolysis  (Reference  32)  and  in  shock  tube 
pyrolysis  (References  33,34).  This  reaction  has  been  known 
throughout  the  20th  Century  and  is  discussed  in  more  detail  in 
Section  IV. 1. a.  The  rate  constant  for  equation  (12)  is  k  -  1.2  x 
1010  exp  (-51520/RT)  sec"1.  Tsang,  fit  al.  (Reference  34)  also 
found  a  third  reaction  attributed  to  isomerization  of  the  C-N02 
bond,  which  is  represented  by  equation  (13)  with  k  estimated  to  be 


1013  exp  (-55980/RT)  sec'1.  However,  a  discrepancy  exists 
concerning  the  relative  importance  of  equations  (11)  and  (12) . 
Equation  (12)  was  found  to  dominate  in  the  studies  of  Tsang,  et 
al.  (References  33,34)  while  equation  (11)  dominated  in  the 
studies  of  McMillen,  Golden,  and  co-workers  (References  32,42). 

In  both  the  laser  pulse  and  the  shock  tube,  the  temperature  of  the 
reactant  is  stepped  in  several  jxsec  to  the  800  -  1000*0  range. 
Ideally,  the  T-jump  has  a  square  waveform,  but  in  reality  the 
shock  front  can  be  inhomogeneous  and  the  laser  beam  does  not  have 
a  flat  energy  cross-section.  The  different  conclusions  about  the 
relative  importance  of  reaction  channels  (11)  and  (12)  was 
proposed  (References  33,34)  to  originate  in  the  uncertainty  of  the 
reaction  temperature.  Although  both  studies  (References  32,34) 
employed  an  internal  standard  to  measure  temperature,  a  100 *c 
higher  temperature  in  the  laser  pyrolysis  studies  could  account 
for  the  different  results. 

Figure  2  shows  Arrhenius  plots  of  equations  (11),  (12)  and 
(13)  using  the  rate  constants  of  Tsang,  fit  al.  (References  33,34) 
for  equations  (12)  and  (13)  and  the  2-NT  data  for  equation  (11). 

Of  the  three  reactions  in  the  gas  phase,  equation  (12)  has  the 
fastest  rate  up  to  950*C.  At  this  temperature  equation  (11) 
becomes  the  fastest  reaction.  A  similar  analysis  of  the  relative 
rates  of  decomposition  channels  can  be  made  of  1,3,5- 
trinitrobenzene  versus  2  4,6-trinitroaniline.  C-N02  homolysis  of 
TNB  equivalent  to  equation  (11)  has  a  rate  constant  (Reference  52) 
of  k  *  1017-2  exp(-67300/RT)  sec"1.  2 , 4 , 6-Trinitroaniline  is  known 
to  decompose  by  two  branches  shown  in  equation  (14)  to  the  furazan 
and  furoxan  derivatives.  These  cyclized  products  are  equivalent 
to  equation  (12)  for  nitrotoluene .  Two  gas  phase  rate  constant 
measurements  are  available  (References  53,84)  for  equation  (14) 
giving  k  -  10®  8  exp(-38500/RT)  sec"1  and  k  =  1011  exp(-41800/RT) 
sec'1.  Figure  3  shows  a  plot  of  these  three  rate  constants.  The 
cyclization  reaction  equation  (14)  is  the  favored  process  below 
475 *C  or  625 *C  depending  on  which  rate  constant  is  chosen  for 
2, 4,6-trinitroaniline.  C-N02  homolysis  is  the  favored  process 
above  these  temperatures. 
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C-N02  HOMOLYSIS 
FURAZAN/FUROXAN  (A) 
FURAZAN/FUROXAN  (B) 
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Figure  3.  Arrhenius  Plots  of  the  Rates  of  Decomposition  of 
MATB  of  in  the  Gas  Phase  by  C-N02  Homolysis  and 
Fura  z  an/ Furoxan  Format ions . 
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The  conclusion  rsachsd  from  the  appearance  of  Figures  2  and  3 
is  that  C-NOj  hoaolysis  is  the  favored  decomposition  channel  of 
nitroarosatic  compounds  in  the  gas  phase  when  the  temperature  is 
high.  The  cyclization  decomposition  channels  dominate  at  lower 
temperature. 

All  of  the  discussion  above  pertains  to  the  gas  phase 
pyrolysis  mechanism.  Reactions  in  the  condensed  phase  can  be 
strongly  affected  by  near-neighbor  interactions.  For  example,  C- 
N0a  homolysis  in  the  condensed  phase  may  be  assisted  by 
intermolecular  H*  atom  transfer  (Reference  48)  to  produce  HONO 
rather  than  N02« .  In  this  case,  a  concerted  process  involving 
both  C-H  and  c-N02  homolysis  might  determine  the  rate  constant. 
Moreover,  during  shock  and  impact  initiation,  the  temperatures 
reached  are  relatively  high.  For  instance,  Wenograd  (Reference  2) 
estimated  that  TNT  and  TNB  reach  1000-1100 *C  during  impact 
initiation.  Consequently,  C-N02  homolysis  is  likely  to  be 
competitive  with  other  decomposition  channels  under  shock  and 
impact  initiation.  In  fact,  a  number  of  studies  conclude  that 
shock  and  impact  initiation  of  nitrobenzene  derivatives  correlate 
with  attributes  of  the  C-N02  bond.  Owens  (Reference  73)  and 
Politzer  and  co-workers  (References  74,75)  correlated  impact 
sensitivity  with  the  electrostatic  potential  at  the  mid-point  of 
the  C-N02  bond.  Delpuech  and  Cherville  (References  117-119) 
correlated  impact  sensitivity  with  an  index  derived  from  the 
polarity  of  the  C-N02  bond  in  the  ground  and  excited  states  of 
nitroaromatic  compounds . 

Assertions  that  the  rate  of  C-N0Z  hoaolysis  is  fundamentally 
important  in  the  impact  and  shock  sensitivity  of  nitroaromatic 
compounds  (References  73,75)  directly  contradict  conclusions  drawn 
from  other  studies.  In  particular,  Kamlet  and  Adolph  (Reference 
24)  were  able  to  separate  the  impact  sensitivities  of 
nitroaromatic  compounds  into  two  categories  —  those  with  or-CH 
bonds  ortho  to  the  -N02  group  on  the  ring  and  those  without  a-CH 
bonds.  These  results  are  discussed  in  greater  detail  in  Section 
IV. 4. a.  They  imply  that  reactions  of  the  a-CH  bond  as  in  equation 
(12)  are  important  during  the  impact  initiation  process.  In  fact, 
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Dacons,  fit  al.  (Reference  20)  proposed  that  bimolecular  oxidative 
insertion  of  O  into  the  CH  bond  was  the  initial  and  RDS.  Bulusu 
and  Autera  (Reference  66)  observed  a  small  positive  deuterium 
isotope  in  the  shock  sensitivity  and  detonation  velocity  when  TNT 
was  compared  to  TNT-o-d3.  Because  a  primary  isotope  effect  is 
also  observed  in  the  low  temperature  decomposition  of  TNT 
(Reference  51)  they  concluded  that  the  RDS  of  slow  decomposition, 
shock  initiation,  and  detonation  are  the  same  and  involves  C-H 
bond  scission.  The  small  size  of  the  isotope  effect  implies, 
however,  that  significant  dilution  has  taken  place  during  shock 
initiation  and  detonation. 

The  direct  conflict  concerning  whether  C-N02  homolysis  or  a- 
CH  bond  scission  is  rate  controlling  in  nitroaromatic  compounds 
may  simply  stem  from  the  temperature  conditions  of  the  experiment. 
As  shown  in  Figure  2,  the  rates  of  C-N02  homolysis  and  the 
cyclization  reaction  equation  (12)  are  competitive  in  the  850  - 
1200*0  range  where  shock  and  impact  initiation  takes  place.  This 
would  imply  that  the  chemical  reactions  responsible  for  shock  and 
impact  initiation  are  not  necessarily  the  same  as  those  of  slow, 
controlled  thermolysis.  In  fact.  Storm,  fi£  al.  (Reference  27) 
found  TNT  to  be  an  outlying  point  in  a  plot  of  impact  vs.  shock 
sensitivity  for  21  explosives. 

In  conclusion,  C-N02  homolysis,  whether  unimolecular  as  in 
the  gas  phase  or  assisted  by  intermolecular  interactions  in  the 
condensed  phase,  is  the  dominant  early  thermal  decomposition 
reaction  of  nitroaromatic  compounds  at  high  temperature,  probably 
above  900 *C. 

2 .  Nitro-Nitrite  Isomerization 

Evidence  of  the  nitro-nitrite  rearrangement  on  an  aromatic 
backbone  exists  primary  in  the  gas  phase  and  in  dilute  solution. 

In  the  condensed  phase  of  the  bulk  material,  the  evidence  is  more 
circumstantial.  Unfortunately,  no  aryl  nitrite  compounds  appear 
to  have  been  isolated  and  characterized. 

In  the  gas  phase  the  clearest  evidence  of  the  nitrobenzene- 
to-nitritobenzene  isomerization  comes  from  the  fact  that  the 
thermal  decomposition  of  nitrobenzene  is  insensitive  to  the 
presence  of  excess  N02  (Reference  32).  An  estimate  (Reference  33) 
of  the  rate  constant  for  this  rearrangement  produces  the  Arrhenius 
plot  shown  in  Figure  2.  While  this  estimate  shows  that 
isomerization  is  never  the  fastest  reaction,  a  small  change  in  the 
assumed  A  factor  alters  the  importance  of  the  isomerization 
process  relative  to  the  other  reactions  in  Figure  2.  The 
isomerization  reaction  also  has  a  small  or  negative  volume  of 
activation,  and,  therefore,  would  be  favored  relative  to 
unimolecular  bond  scission  at  high  pressure  (Reference  120) . 
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In  the  condensed  phase,  isomerization  would  be  expected  to 
increase  in  importance  relative  to  the  gas  phase  because  of  the 
low  volume  of  activation.  Evidence  of  a  reversible  nitro-nitrite 
rearrangement  involving  a  nitroaromatic  compound  exists  in  the  *H 
and  UC  NMR  spectra  of  2 , 6-dichloro-4-methylphenol  following 
nitration  in  chloroform  (Reference  121) .  Nitration  occurs  at  the 
4-position  which  makes  this  carbon  atom  formally  sp3  hybridized. 
The  nitro-nitrite  rearrangement,  therefore,  does  not  actually 
occur  on  an  aromatic  carbon  atom.  Nitro-nitrite  isomerization  has 
long  been  known  to  take  place  on  an  aliphatic  carbon  atom 
(Reference  122). 

A  strong  item  of  circumstantial  evidence  for  the  nitro- 
nitrite  rearrangement  of  nitroaromatic  compounds  in  the  bulk  state 
is  the  fact  that  N02(g)  is  not  observed  or  is,  at  most,  a  minor 
product  of  thermal  decomposition  below  500 *C  (Reference  123) . 
Instead,  a  large  quantity  of  NO  is  formed  (References  48,49,123- 
125) .  These  observations  have  been  cited  as  evidence  of  the 
different  contributions  of  equations  (15)  and  (16)  (Reference 
125);  (17)  (References  32,33,49,120);  and  (18)  (References 
33,49,125,126) . 


As  with  equation  (15)  in  Section  IV. 1.,  equations  (16) -(18)  have 
been  analyzed  most  extensively  in  the  gas  phase,  but  not  in  the 
condensed  phase,  where  the  details  are  more  difficult  to 
establish.  Beynon,  al .  (Reference  127)  were  among  the  first  to 
contribute  to  this  area.  They  observed  that  pyrolysis  of 
nitrobenzene  at  600 *C  produced  NO  and  the  phenoxy  radical.  Above 
800 *C,  shock  tube  studies  (Reference  33)  indicated  that  equations 
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(17)  and  (18)  occurred,  while  laser  powered  pyrolysis  revealed 
(Reference  32)  that  the  nitrite  fores  by  equation  (17),  rather 
than  equations  (15)  and  (16).  However,  equation  (15)  is  more  than 
twice  as  prevalent  as  equation  (17)  above  800*C  (Reference  33). 

The  fact  that  these  shock  tube  and  laser  pyrolysis  studies  were 
conducted  at  low  pressure  would  disfavor  the  bisolecular  equation 
(16).  Equation  (16)  would  be  expected  to  be  such  sore  important 
in  the  condensed  phase. 

Values  of  Ea  and  the  BDE  for  equation  (15)  are  in  the  65-72 
kcal/mol  range  (Section  IV. 1.).  Although  the  rate  of 
isomerization  by  equation  (17)  has  not  been  measured  directly,  it 
is  likely  to  be  slow  based  on  the  estimate  of  a  relatively  high  E, 
of  55.9  kcal/mol  (Reference  33)  and  computational  estimates 
(Reference  120).  The  BDE  of  equation  (18)  is  approximately  23.7 
kcal/mol  (Reference  128) .  Taken  together  these  results  are 
consistent  with  equation  (15)  being  the  prevalent  reaction  at  high 
temperature  and  low  pressure  (Reference  33) . 

Reactions  of  the  type  (15)  -  (18)  also  appear  to  occur  when 
other  substituents  are  present  on  the  aromatic  ring.  Field  and 
Meyerson  (Reference  116)  found  ions  in  the  mass  spectrum 
consistent  with  equation  (18)  when  1,3-  and  1,4-nitrotoluene  and 
1,3-dinitrobenzene  decompose.  The  major  pathway  of  decomposition 
of  3,4,5-trinitrotoluene  in  the  EI/MS  is  loss  of  NO,  which  was 
proposed  to  result  from  equations  (17)  and  (18)  (Reference  39). 
Decanitrobiphenyl  was  speculated  to  decompose  in  solution  through 
the  isomerization  steps,  equation  (17)  (Reference  129).  MNDO  and 
MINDO/3  computations  (Reference  130)  on  a  model  compound  of  TNT, 
1-nitropropene,  reveal  that  the  transition  state  equivalent  to 
equation  (17)  is  favorable  and  is  consistent  with  the  formation  of 
radicals  by  equation  (18) . 

Considering  the  evidence  available,  concrete  proof  does  not 
exist  for  nitro-nitrite  isomerization  of  nitroaromatic  compounds 
during  thermal  decomposition.  However,  circumstantial  evidence  is 
abundant.  The  isomerization  equation  (17)  could  occur  directly  in 
the  condensed  phase,  or  at  higher  temperatures  it  could  occur  by 
equations  (15)  and  (16).  A  variation  of  equation  (16),  for  which 
there  is  abundant  evidence  (References  121,122,131),  involves 
converting  an  sp2  carbon  to  an  sp3  carbon  and  adding  N02  at  the  sp3 
carbon  atom.  The  copious  amount  of  N0(g)  formed  upon  thermolysis 
of  nitroaromatic  compounds  is  further  evidence  of  equation  (18). 
However,  it  is  doubtful  that  all  of  the  NO(g)  generated  by 
nitroaromatic  compounds  comes  from  equation  (18) .  Other  sources 
of  NO(g)  need  to  be  considered.  In  addition  to  equation  (18), 
N0(g)  can  be  formed  by  thermal  decomposition  of  HONO,  which  forms 
when  •  N02  and  H*  adventitiously  encounter  one  another.  Most 
nitroaromatic  compounds  are  abundant  sources  of  *H  (Reference  61). 
Equation  (19)  occurs  readily  at  the  temperature  which  nitrobenzene 
compounds  thermally  decompose.  This  route  to  N0(g)  is  indicated 
to  be  unimportant  during  pyrolysis  of  nitrobenzene  in  one  study. 
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2HQNO  —  NO  *  NC^  +  Kfi 


(It) 


(Reference  49)  but  contributing  in  another  (Reference  125) . 
Equation  (19)  has  been  mentioned  frequently  in  the  thermal 
decomposition  mechanisms  of  substituted  nitroaromatics  (References 
33,36/48,130,132) .  A  third  source  of  NO(g)  is  thermal 
decomposition  of  nitrosoaromatic  compounds  that  form  in  the 
decomposition  scheme  of  nitroaromatic  derivatives  (Section  IV. 3). 

3.  Formation  of  Nitrosoaromatic  intermediates 

Kinetics  measurements  of  the  formation  and  thermal 
decomposition  of  nitrosoaromatic  compounds  from  nitroaromatic 
compounds  are  not  available.  However,  from  a  thermochemical  point 
of  view,  the  C-N  BDE  of  nitrosobenzene  is  50.811  kcal/mol  compared 
to  nitrobenzene  at  71.311  kcal/mol  (Reference  133).  On  this 
basis,  formation  and  reaction  of  nitrosoaromatic  compounds  is 
probably  a  lower  temperature  process  than  C-N02  homolysis.  Any 
temperature  at  which  C-N02  homolysis  is  facile  exceeds  the 
temperature  required  to  cleave  the  C-NO  bond. 

As  is  the  case  for  other  reactions  of  the  -N02  substituent 
described  in  Sections  IV. 1.  and  2.,  the  clearest  evidence  of 
formation  of  the  nitroso  substituent  from  the  nitro  group  exists 
in  the  gas  phase.  Equations  (20)  -  (23)  give  several  possible 
routes  to  the  nitroso  group. 
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Equation  (20)  has  boon  proposed  on  the  basis  of  pyrolysis  of 
nitrobenzene  in  benzene  solution  at  275 *C  where ,  after 
chromatography ,  nitrosobenzene  is  the  dominant  nitrogen-containing 
product  (Reference  49) .  In  the  EI/KS  (Reference  39)  and  MS /MS- 
CZD  (Reference  134)  of  TNT,  an  ion  with  m/e  -211  corresponding  to 
the  loss  of  0  is  important. 

Equation  (21) ,  in  which  an  oxygen  atom  is  abstracted  by 
reaction  with  a  radical,  would  be  expected  with  a  radical- 
producing  substituent  ortho  to  the  -N0a  group.  Since  it  can  be  a 
bimolecular  reaction,  equation  (21)  is  potentially  important  in 
the  condensed  phase.  Equation  (21)  has  been  proposed  to  occur  for 
nitrobenzene  in  benzene  solution  (Reference  49)  and  in  the  gas 
phase  of  TNT  by  EI/MS  (Reference  39),  and  MS/ MS -Cl D  (Reference 
134)  where  loss  of  OH  produces  the  dominant  ion. 

The  fact  that  NO  forms  in  high  concentration  when  many  bulk 
nitroaromatic  compounds  are  pyrolyzed  (References  48,49,123-125) 
makes  equation  (22)  a  probable  source  of  nitrosoaromatic  compounds 
by  recombination  in  the  condensed  phase.  Equation  (23)  has  been 
analyzed  computationally  as  a  potential  step  in  the  decomposition 
of  TNT  (Reference  135) .  MINDO  calculations  suggest  that  equation 
(23)  is  favored  over  H  transfer  from  -CH3  to  -N02  to  form  the  aci- 
isomer  (Reference  135) .  The  same  conclusion  was  reached  from  MNDO 
and  MINDO/ 3  calculations  of  1-nitropropene  (Reference  130)  a  model 
compound  for  TNT.  Oxidative  insertion  of  0  into  the  CH  bond  has 
been  speculated  upon  as  the  RDS  in  the  decomposition  of  TNT 
(Reference  20) . 

The  chemistry  of  nitrosoaromatic  derivatives  must  be 
considered  in  the  lower  temperature  decomposition  schemes  of 
nitroaromatic  explosives,  but  probably  not  in  shock  initiation 
processes.  Evidence  for  this  comes  from  XPS  studies  of  Owens  and 
Sharma  (Reference  136)  who  found  evidence  of  the  nitroso  group  in 
thermally  decomposed  TNT,  but  not  in  shocked  TNT.  Indeed,  Adams, 
et  al.  (Reference  137)  isolated  a  toluene-soluble  nitroso- 
containing  compound  proposed  to  be  2-nitroso-4 , 6-dinitrobenzyl 
alcohol  from  the  thermal  decomposition  products  of  TNT.  Thermal 
decomposition  of  TNT  alone  and  in  the  presence  of 
hexamethylbenzene  produces  a  nitroxide  radical  that  corresponds  to 
condensation  with  loss  of  OH,  equation  (24)  (Reference  138) .  This 
radical  is  closely  related  to  the  nitroso  compounds  discussed 
above.  Likewise,  nitrosoaromatics  can  couple  to  form  azoxy  and 
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•so  compounds  which  have  been  isolated  fro*  th«  residue  of  TNT 
for— d  during  Id  hrs  at  200*C  (Reference  20) .  Coupling  of 
nitrosoaroaatics  to  fora  aso  and  azoxy  compounds  is  proposed  in 
ottear  studies  (Reference  139) .  Equation  (25)  is  representative. 

Decomposition  of  nitrosoaroaatic  compounds  is  a  source  of  NO 
upon  thermolysis.  Both  nitrosobenzsne  (Reference  139)  and  the 
nitroso  products  of  TNT  (Reference  134)  produce  NO  by  MS/KS-CID. 


4.  Reactions  of  the  -N02  Group  with  an  ortho  Substituent 

Numerous  studies  indicate  that  most  polyatomic  substituents 
positioned  ortho  to  the  -N02  group  on  the  aromatic  ring  enhance 
the  reactivity  compared  to  mate  and  para  positioning.  Synthesis 
chemists  knew  this  fact  in  the  19th  Century.  Reactions  of 
substituents  with  an  ortho  -N02  group  are  now  a  widely  used 
synthetic  route  to  many  heterobicyclic  molecules  (References  140- 
142) .  Mechanisms  of  these  reactions  have  been  discussed  in  the 
broad  sense  (References  140,143).  Reactivity  differences  in 
relation  to  explosive  properties  began  to  be  discussed  in  the 
early  20th  Century  (References  83,144).  For  example,  2,4,6-TNT 
self-ignites  at  about  210 *C,  whereas  1,3,5-TNB  is  stable  at  300 *C 
indicating  that  -CH3  somehow  activates  the  molecule. 

Clear  illustrations  of  the  generality  of  the  "ortho  effect" 
arise  when  classes  of  substituted  nitrobenzene  compounds  are 
compared.  The  relative  rates  of  thermal  decomposition  of  2- 
nitrotoluene,  3-nitrotoluene,  and  4-nitrotoluene  in  the  vapor 
phase  at  200*C  are  400:0.8:2.0  (Reference  96).  The  EI/MS  and 
MS/MS-CID  of  these  molecules  differ  significantly.  2 -NT  loses  OH 
while  3 -NT  and  4 -NT  lose  0,  NO,  and  N02  (References  36,126). 

-CH3* • • ON-  Intermolecular  interaction  clearly  destabilizes  the 
molecule. 

Kinetic  data  for  the  gas  phase  decomposition  of  a  series  of 
ortho-X-nitrobenzene  derivatives  (Reference  102)  in  Table  8  are 
further  evidence  of  the  generality  of  the  ortho  effect.  In  all 
three  compounds  the  primary  deuterium  isotope  effect  indicates 
that  the  H  atom  is  transferred  in  the  transition  state  of  the 
reaction.  Ea  for  the  2 -NT  in  Table  8  is  10  -  20  keal/mol  lower 
than  those  in  Table  7  where  the  C-N02  homolysis  occurs.  All  three 
compounds  in  Table  8  primarily  fragment  by  loss  of  OH  in  the  EI/MS 
and  MS/MS-CID  modes  (References  36,126). 


The  reactions  between  -NO,  and  ortho  related  substituents 
profoundly  affect  the  thermal  decomposition  chemistry  and 
explosive  properties  of  polynitroaromatic  explosives.  Individual 
substituents  are  now  considered  by  category. 

a.  a-CH  Bonds 

The  number  of  clear  indications  are  overwhelming  that  an  a-CH 
bond  ortho  to  -NO,  on  an  aromatic  ring  activates  the  thermal 
decomposition  process  (References  5,20,21,24,33,34,39,51,53,86, 
96,126,132,134,140,145-147).  That  an  alkyl  group  plays  a 
fundamental  role  in  thermal  decomposition  is  not  at  all  surprising 
when  these  compounds  are  viewed  as  a  subclass  of  hydrocarbon 
chemistry.  The  RDS  in  the  oxidation  of  hydrocarbons  in  the  liquid 
phase  involves  C-H  scission  (Reference  148).  The  C-H  bond  of  an 
alkyl  group  is  more  reactive  than  a  C-H  bond  of  an  aryl  group 
(Reference  148) .  However,  for  destabilization  of  the  nitrobenzene 
fragment,  it  is  not  simply  the  presence  of  an  ortho  alkyl  group 
that  counts.  The  presence  of  an  a-CH  bond  in  the  ortho  alkyl 
group  is  the  most  significant  factor.  Evidence  for  this  is  the 
fact  that  at  200  -  260  *C  in  the  melt  phase,  2,4,6-TNB  substituted 
by  -CH3,  -CH,CH3 ,  and  -CH(CH3) ,  decompose  at  about  the  same  rate, 
whereas  -C(CH3)3  substitution  produces  a  30  -  50  times  slower  rate 
(Reference  86) .  Following  a-CH  abstraction  the  reaction  sequence 
for  this  series  of  compounds  is  very  different  (Reference  37). 

Chemical  details  about  how  a  nitroaromatic  compound  becomes 
activated  by  ortho  -CH3  substitution  cause  to  light  nearly  a 
century  ago  when  2-NT  was  found  to  condense  to  anthranil,  equation 
(12) ,  upon  heating  (References  149-151) .  In  particular,  2-NT 
heated  in  aqueous  or  alcoholic  base  rearranged  to  anthranilic  acid 
(References  149-151)  which  is  the  hydrolysis  product  of  anthranil. 
Likewise,  Fields  and  Neyerson  (Reference  116)  recovered 
methylanthranilate  when  2-NT  in  CH30H  was  heated  to  600  *C  in  a 
vycor  tube  for  11  sec.  The  mechanism  of  this  reaction  has  been 
discussed  (Reference  65) .  Because  these  reactions  occur  in  the 
condensed  phase,  the  role  played  by  bimolecular  processes  and  wall 
effects  is  unknown.  These  perturbations  are  eliminated  by  the  use 
of  a  single  pulse  shock  tube  (References  33,34).  The  kinetics  and 
products  of  unimolecular  thermal  decomposition  of  2-NT  reveal  that 
equation  (12)  is  the  dominant  reaction  pathway  at  least  below 
about  900*0  (References  33,34).  Previous  kinetic  studies  of  2-NT 
using  the  rate  of  gas  evolution  (References  96,102)  and  these 
shock  tube  data  are  given  in  Table  9.  The  Arrhenius  factors  are 
smaller  than  those  for  C-N02  homolysis  (Table  7) .  The  relatively 
low  A  factor  suggests  a  tight  transition  state  that  is  expected  of 
a  cyclic  activated  complex  collapsing  to  a  closed  ring  (References 
34,84).  Indeed,  AS*  is  -15  cal/deg  for  2-NT  compared  to  20  and 
14.5  cal/deg  for  3 -NT  and  4 -NT  (Reference  96).  As  described  in 
Section  IV. 1,  3 -NT  and  4 -NT  thermally  decompose  by  C-N02  homolysis 
which  produces  a  positive  AS*.  The  clear  experimental  evidence 
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Table  8.  Gas  Phasa  Kinatic  Data  (Reference  102)  for 
Nitrobenzene  Darivativas  having  ortho-X 
Substitution 


X 

T,  *C 

Ea,  kcal/sol 

CHS 

350  -  420 

1.54 

49.5 

HHa 

400  -  470 

1.48 

57.5 

OH 

390  -  460 

1.75 

52.4 

Tabla  9.  Arrhanius  Data  for  Condensation  of  2-Nitrotoluene 
to  Anthranil  (aquation  12) 


T,  *C 

Ea,  kcal/nol 

A,  sec*1 

Ref. 

300  -  350 

42.5 

lO10.22 

96 

400  -  450 

49.5  ±  1.3 

1012.*  ±  0  * 

102 

797  -  907 

51.520 

1013.08 

33 
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that  2-NT  condansas  to  anthranil  upon  heating  is  support ad  by 
tharalativa  anargiaa  determined  by  ab  initio  SCF-MO  methods 
(Reference  152) .  The  bicyclic  tautoaar  intemediate  and  the 
anthranil  product  of  aquation  (26)  fora  exotheraically  and  thus 
are  theraodynaaically 


favored  relative  to  2-NT.  The  first  step  of  equation  (26)  is 
suggested  by  UV  photolysis  of  2-NT  which  produces  a  short-lived 


(D 


aci  tautoaer  of  2-NT  (I) ,  (References  153-156) .  (I)  is  the 

precursor  of  the  intemediate  species  in  equation  (26) .  SCF-MO 
calculations  on  2-NT  at  the  STO-3G  level  indicate  that  the  barrier 
to  H  transfer  when  foraing  (I)  is  54  kcal/mol  (Reference  157) , 
which  is  4-13  kcal/mol  larger  than  the  global  induction  phase 
activation  energy  of  2-NT  (Table  5) .  An  earlier  MNDO  calculation 
on  TNT  concluded  (Reference  135)  that  intranolecular  H  transfer  to 
font  (I)  was  less  likely  than  0  insertion  into  the  CH  bond  by 
equation  (23) .  However,  the  SCF-MO  aodel  is  nore  dependable  for 
assessing  the  relative  energy  of  these  processes.  Elsewhere  the 
MNDO  model  of  1-nitropropene  indicated  that  H  transfer  from  -CH3 
to  -N02  is  energetically  favorable  (Reference  130) .  This  result 
is  consistent  with  the  SCF-MO  results. 

The  experimental  and  computational  results  for  2-NT  are  useful 
guides  to  the  possible  behavior  of  TNT.  However,  2-NT  is  only  a 
model  of  an  explosive  and  has  been  studied  primarily  in  the  gas 
phase.  From  the  point  of  view  of  an  explosive,  the  behavior  of 
the  bulk  material  in  the  condensed  phase  must  be  understood. 
Consequently,  the  unimolecular  reaction  conditions  of  the  gas 
phase  need  to  be  modified  to  account  for  bimolecular  reactions  and 
catalysis  that  are  important  in  the  condensed  phase.  Condensed 
phase  decomposition  is  further  complicated  by  the  fact  that  H 
atoms  are  such  important  participants  in  the  decomposition  of  TNT 
(References  61,158).  The  rate  of  H  transfer  is  fast,  being  less 
than  3.1  psec  in  solution  (Reference  156).  Consequently,  many  H 
atom  transfer  steps  can  occur  in  a  very  short  time  at  temperatures 
below  300 *C  where  most  decomposition  studies  of  polynitroaromatic 
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explosives  are  conducted.  Moreover,  0  transfer  occurs  under 
certain  conditions  in  the  condensed  phase  (Reference  159). 

TNT  selts  at  81 *C  so  that  studies  of  condensed  phase 
decomposition  almost  always  refer  to  the  molten  state.  Several 
studies  report  separation  and  identification  of  products  in  an 
attempt  to  determine  how  the  methyl  group  activates  thermal 
decomposition.  Many  products  are  possible  because  all  of  the 
oxidized  forms  of  -CH3  and  reduced  forms  of  -N02  are  present  and 
each  displays  its  own  intermolecular  and  intramolecular  chemistry. 
Consequently,  anthranil,  alcohol,  aldehyde,  carboxylic  acid, 
nitroso,  oxime,  nitrile,  nitrone,  nitroxide,  azo,  and  azoxy  groups 
form  (References  20,21,137,139,160-162).  Disruption  of  the 
aromatic  ring  also  occurs  leading  to  gas  products  including  NO, 

CO,  C02,  N2  and  H20  (References  87,137).  An  insoluble  high 
molecular  weight  "explosive  coke"  having  the  elemental  composition 
of  C«H3N203.7S  forms  from  TNT  (Reference  137) .  All  of  these 
products  result  from  "cooking"  reactions  and,  as  such,  can  not  be 
attributed  with  confidence  to  the  initiation  or  explosion 
chemistry  of  TNT.  However,  several  useful  results  emerge  from 
this  work  and  are  discussed  below. 

Rogers  (Reference  21)  used  thin-layer  chromatography  to 
separate  the  products  formed  by  programmed  heating  of  TNT  at  300 
*C.  1,3,5-TNB,  2,4-dinitroanthranil,  2 , 4 , 6-trinitrobenzyl 

alcohol,  and  2 , 4 , 6-trinitrobenzoic  acid  were  isolated  from  a  large 
quantity  of  unreacted  TNT.  He  proposed  that  oxidation  of  -CH3 
to  -CHjOH  is  the  first  step.  Dacons,  fifc  al.  (Reference  20)  heated 
TNT  for  16  hrs  at  200 *c  in  a  loosely  capped  vial.  About  25 
benzene-soluble  compounds  were  extracted  from  the  dark  brown 
residue.  An  insoluble  residue  included  the  "explosive  coke" 
reported  by  Adams,  fit  al.  (Reference  137).  Unreacted  TNT 
constituted  75  -  90%  of  the  material.  The  same  benzene-soluble 
products  found  by  Rogers  (Reference  21)  (except  for  TNB)  were 
obtained,  but,  in  addition,  2 , 4 , 6-trinitrobenzaldehyde  and  dimeric 
reduction  products  suspected  of  containing  azo  and  azoxy  linkages 
(II)  were  found.  None  of  the  products  amounted  to  more  than  1%  of 


(H) 


the  original  quantity  of  TNT.  Neither  Rogers  (Reference  21)  nor 
Dacons,  fit  al.  (Reference  20)  analyzed  the  gases  formed,  but 
Rogers  noted  that  only  a  few  percent  of  the  weight  is  lost. 
Consequently,  most  of  the  decomposition  reactions  produced 
products  with  low  volatility.  It  is  not  surprising  that  the 
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products  in  these  two  studies  differ  somewhat  because  of  the 
different  reaction  temperature.  While  Dacons,  fit  al.  (Reference 
20)  stress  the  importance  of  these  differences,  Rauch,  fit  ala 
(References  163-165)  repeated  this  work  by  using  TNT  heated  at 
275*C,  which  is  intermediate  between  the  temperatures  used  by 
Rogers  (Reference  21)  and  Dacons,  fit  al.  (Reference  20) .  Products 
common  to  both  studies  were  isolated  indicating  that  the  main 
factor  in  the  product  differences  was  simply  the  reaction 
temperature.  After  merging  the  results  of  these  cooking 
degradation  studies,  the  presence  of  successive  oxidation  products 
of  the  -CH3  group,  equation  (27) ,  is  conspicuous.  All  of  these 
products  also  appear  in  the  synthesis  of  TNT. 

Many  of  the  products  discussed  above  are  formed  from  TNT 
because  of  the  relatively  low  decomposition  temperature  and  the 
long  reaction  time.  It  is  unlikely  that  they  all  are  present  or 
contribute  significantly  during  rapid  decomposition  where  the 
reaction  times  are  shorter  and  the  temperature  is  higher.  As  has 
already  been  emphasized,  high  temperature  favors  higher  Ea  and 
more  endothermic  processes.  Therefore,  experiments  designed  to 


specify  individual  steps  or  to  control  the  conditions  more 
precisely  require  interpretative  care  if  the  goal  is  to  define  the 
decomposition  mechanism  of  TNT.  For  example,  a  low  Ea  process  is 
revealed  by  ESR  spectroscopy  during  UV  photolysis  of  TNB  in 
toluene  at  200  K.  The  primary  radical  is  (III)  which  converts 
with  time  or  on  heating  to  two  conformers  of  (IV)  (Reference  166) . 
Caution  must  be  exercised  when  using  these  results  to  understand 
TNT  decomposition  because  these  are  low  Ea  processes.  Moreover, 
Davies,  fit  al.  (Reference  167)  found  that  the  radical  formed  by 
photolyzing  TNT  in  diethylether  is  not  the  same  as  those  formed  by 
heating  TNT.  On  the  other  hand,  there  is  absolutely  no  doubt 
about  the  ease  by  which  an  H  atom  can  migrate  from  the  -CH3  group 
of  TNT.  In  solution,  UV  photolysis  gives  evidence  of  H  atom 
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(HI)  (IV) 

transfer  to  give  (I')  (References  153-156),  which  is  in 
equilibrium  with  the  trinitrobenzyl  anion  (V)  (equation  28) 
(References  145,156).  Lewis  bases,  such  as  OCH3‘  and  OC2H5' , 


O’)  oo 

abstract  IT  from  the  -CH3  group  of  TNT  (References  158,168)  to 
produce  (V) .  (V)  is  highly  reactive  and  under  appropriate 

conditions  is  able  to  couple  to  produce  HNS  (Reference  169) .  In 
the  gas  phase,  the  dominant  fragmentation  step  of  TNT  and  1,3,5- 
trinitroxylene  is  loss  of  OH  (References  39,126,132,134,170).  For 
this  to  occur  H  atom  transfer  to  -N02  followed  by  -N(0)-0H 
homolysis  must  be  the  dominant  initial  decomposition  step  of  the 
isolated  molecule.  In  fact,  all  TNT  isomers,  except  for  3,4,5- 
TNT,  fragment  by  loss  of  OH  (Reference  146) . 

Solid  evidence  of  C-H  bond  scission  from  the  -CH3  group  of 
neat  TNT  is  the  primary  deuterium  isotope  effect  observed  by  DSC 
during  the  induction  phase  (Figure  1)  of  TNT  and  TNT-a-d3 
(Reference  51) .  A  primary  isotope  effect  indicates  that  C-H  bond 
scission  occurs  in  the  RDS.  Shackelford,  al.  (Reference  51) 
proposed  that  a  single  unidentified  catalyst  promotes  the  reaction 
rate  in  the  induction  phase.  They  were  apparently  unaware  of  the 
earlier  report  of  Maksimov,  al.  (Reference  59)  who  suggested 
that  2,4-DNAn  or  2 , 4 , 6-trinitrobenzaldehyde  (equation  (27))  was 
the  catalyst.  Maksimov,  fijfc  al. .  found  that  a  9:1  TNT: 2,4-DNAn 
mixture  reproduced  the  acceleratory  decomposition  characteristics 
of  TNT  (Reference  59).  2 , 4 , 6-Trinitrobenzaldehyde  was  shown 

previously  to  catalyze  the  decomposition  of  TNT  (Reference  20) . 

In  reality  it  is  very  difficult  to  specify  the  catalysts  of 
decomposition  of  TNT  with  certainty.  Robertson  found  that  all  25 
compounds  that  he  added  to  TNT  were  catalysts  (Reference  44) . 
Molecules  that  are  able  to  produce  H*  catalyze  the  decomposition 
(Reference  51) .  Unfortunately,  ample  concentrations  of  such 
species  (10~*  M)  appear  at  relatively  low  temperature  when  neat 
TNT  is  heated  (Reference  61) . 
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Whether  it  is  intra-  or  intermolecular  H*  transfer  that 
produces  the  aci  tautomer  of  TNT  (I1)  and  starts  the  decomposition 
chain  is  more  of  academic  mechanistic  interest  than  practical 
importance.  However,  some  studies  emphasize  intra-  or 
intermolecular  H-  transfer  (Reference  48)  to  the  exclusion  of  the 
other.  Considering  all  of  the  literature  together,  there  is 
little  doubt  that  when  both  are  possible,  both  will  contribute. 

In  solution  and  in  the  gas  phase  at  low  pressure,  unimolecular  H* 
transfer  dominates  (References  32,34,36,37,96,102,116,155,156). 

If  the  pressure  is  raised  in  the  vapor  phase,  the  decomposition  of 
TNT  shifts  from  unimolecular  at  lower  pressure  to  bimolecular  at 
higher  pressure  (Reference  53).  In  the  face  of  this  result  it  is 
interesting  that  the  decomposition  rate  of  molten  TNT  is 
independent  of  pressure  up  to  50  kbar  (Reference  171) .  This  has 
been  attributed  to  the  fact  that  the  volume  of  activation  of  2,4- 
DNA  is  similar  to  that  of  TNT  (Reference  172) .  On  the  other  hand, 
it  could  mean  that  intermolecular  association  and  H*  transfer 
dominate  at  all  pressures  studied. 

Evidence  of  significant  intermolecular  association  of 
polynitrobenzene  molecules  in  the  condensed  phase  that  would 
facilitate  intermolecular  H*  transfer  appears  in  many  forms.  The 
crystal  structure  of  TNT  (References  173,174)  reveals  that 
intramolecular  and  intermolecular  H*  transfer  may  be  equally 
probable  in  the  condensed  phase.  The  intermolecular  -CH* • • ON- 
distances  (2.39  -  2.62A)  are  very  similar  to  the  intramolecular 
distances  (2.38  -  2.70A).  Structural  details  of  the  solid  phase 
do  not  necessarily  carry  over  to  other  phases,  but  the  vibrational 
spectra  of  TNT  in  the  solid,  liquid,  and  vapor  state  are  similar 
(References  174-177) .  This  suggests  that  the  intermolecular 
-CH* • • ON-  distances  in  the  solid  are  similar  to  those  in  the  melt 
where  decomposition  occurs.  The  large  amount  of  intermolecular 
association  in  the  melt  is  responsible  for  the  unusually  high 
Trouton's  ratio  of  AHvap/Tb  =  23  -  33  cal/deg*mol  (Reference  48) 
for  polynitrobenzene  compounds  compared  to  about  22  cal/ deg* mol 
for  weakly  associating  molecules.  Aggregates  of  nitrobenzene 
molecules  exist  even  at  350 *C,  (Reference  178).  The  A  factors  for 
thermal  decomposition  of  many  polynitroaromatic  compounds 
(Reference  48)  are  in  the  range  of  105  -1010  sec'1  which  is  smaller 
than  the  lower  limit  of  unimolecular  reactions  of  about  1012  sec'1. 
The  negative  AS*  for  some  of  these  reactions  is  suggestive  of  an 
ordered  transition  state  (Reference  34) .  An  energetically  stable 
transition  state  forms  by  intermolecular  H*  transfer  from  -CH3 
to  -no2  according  to  MNDO  and  MINDO/3  MO  calculations  on  1- 
nitropropene  (Reference  130) .  Finally,  the  mere  fact  that  the 
rate  of  thermal  decomposition  of  TNT  is  about  ten  times  faster  in 
the  melt  than  in  the  vapor  phase  (Reference  48)  is  reasonable 
evidence  of  the  importance  of  intermolecular  H*  transfer  in  the 
condensed  phase.  However,  this  latter  observation  could  also  be 
attributed  to  catalysis  by  the  decomposition  products. 
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C-H  bond  scission  from  the  -CH3  group  of  TNT  appears  to  play  a 
practical  role  in  the  sensitivity  and  explosive  properties  of 
polynitroaromatic  explosives  in  the  bulk  state.  The  explosive 
properties  of  interest  include  the  sensitivity  to  impact  and  shock 
initiation.  As  discussed  in  Section  II. 1,  impact  sensitivity  is 
determined  by  the  drop  weight  impact  test  and  is  governed  to  a 
large  extent  by  thermal  decomposition  processes  (References  1,2). 
Thermal  decomposition  is  influenced  by  the  structure  and 
composition  of  the  molecule.  In  accordance  with  this  notion,  log 
Hso  for  polynitroaromatic  compounds  is  found  to  correlate  with  the 
oxidant  balance  (Reference  24) ,  the  reducing  and  oxidizing 
valencies  of  the  composing  elements  (Reference  179) ,  the  molecular 
electronegativity  (Reference  180) ,  and  the  polarity  of  the  C-N02 
bond  in  the  ground  and  excited  states  (Reference  117) . 

Particularly  interesting  results  are  those  of  Kamlet  and  Adolph 
(Reference  24)  and  Bliss,  al.  (Reference  181)  who  found  a 
statistically  significant  separation  of  log  H50  vs.  OB100  for 
nitroaromatic  explosives  into  two  families.  With  several 
exceptions  the  distinguishing  feature  is  the  presence  or  absence 
of  an  cr-CH  bond  ortho  to  -N02.  Compounds  with  a-CH  bonds  were 
more  sensitive  to  impact.  Bliss,  al.  (Reference  181)  noted 
that  this  separation  may  merely  be  the  result  of  dilution  of  OB100 
by  the  CH3  group.  Bulusu  and  Autera  (Reference  66)  noted  a  small 
deuterium  isotope  effect  in  the  shock  sensitivity  (1.09  -  1.14) 
and  detonation  velocity  (1.04  -  1.06)  when  TNT  was  compared  to 
TNT-a-d3 .  This  observation  suggests  that  a-CH  bond  scission  might 
have  some  role  in  shock  initiated  chemistry  and  in  steady 
detonation.  Consistent  with  this  is  the  fact  that  the  shifts  of 
the  N(ls)  and  O(ls)  XPS  peaks  in  shocked  TNT  have  been  interpreted 
as  showing  migration  of  H  toward  O  (Reference  136) . 

All  of  the  observations  in  the  foregoing  paragraph  are 
phenomenological  in  that  they  suggest  a  role  for  the  a-CH  bond  in 
the  explosive  properties,  but  do  not  give  a  mechanism.  In  some 
cases  different  processes  produce  the  same  effect.  The  primary 
kinetic  isotopic  effect  (kg/kp  =  1.66  ±  0.2)  found  for  the 
decomposition  of  TNT  over  minutes  at  263 ’C  is  likely  to  result 
from  the  formation  of  2,4-DNA  or  a  related  compound  (Reference 
51) .  The  small  values  for  the  deuterium  isotope  effect  observed 
for  shock  initiation  and  the  detonation  velocity  (Reference  66) 
could  result  from  dilution  of  the  C-H  homolysis  isotope  effect  by 
C-N02  homolysis.  Shock  initiation  and  the  velocity  of  detonation 
are  dominated  by  reactions  in  the  submicrosecond  and  the  higher 
temperature  regime.  In  this  regime  an  isotope  effect  can  also 
arise  from  the  transport  rate  differences  of  H  and  D  (1:1.4). 
Therefore,  the  small  isotope  effects  observed  may  be  caused  by 
transport  rate  differences  mixed  with  C-H  homolysis  chemistry  and 
diluted  by  C-NOz  homolysis.  Furthermore,  polynitroaromatic 
compounds  containing  a-CH  bonds  are  not  distinguishable  from  other 
compounds  when  the  detonation  pressure  and  Ea  of  slow  thermal 
decomposition  are  correlated  (Reference  7) .  In  fact,  Kamlet  and 
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Hurvitz  (Reference  182)  do  not  believe  that  reaction  kinetics  even 
influences  the  velocity  of  detonation. 


The  overall  description  of  the  decomposition  process  of  neat 
TNT  in  relation  to  the  explosive  characteristics  requires  both 
kinetic  and  mechanistic  footing.  The  following  view  of  TNT  rests 
on  firm  ground  based  on  the  foregoing  discussion: 

(a)  The  dominant  initiating  processes  of  TNT  are: _ -CH, 

oxidation  producing  2.4-dinitroanthranil.  2.4,6- 
trinitrobenzaldehvde .  etc.:  C-NO,  homolysis;  and  catalyst- 
dominated  chemistry.  Validation  of  this  statement  emerges  from 
consideration  of  the  kinetics  and  the  mechanistic  discussions  in 
this  section.  The  global  reaction  rates  for  the  induction, 
acceleratory,  and  decay  phases  of  TNT  below  400*C  were  given  in 
Table  5,  and  for  nitrobenzene  and  nitrotoluene  derivatives  in 
Tables  6,  8  and  9.  With  these  data  the  reaction  time  at  a  given 
temperature  can  be  assessed  for  these  main  classes  of  reactions. 

Ea  values  in  Table  5  for  the  induction  phase  of  TNT  are  41-51 
kcal/mol,  whereas  for  anthranil  formation  from  2-NT  in  Table  9,  E. 
is  the  range  of  42.5  -  51.5  kcal/mol.  This  overlap  is  reasonably 
good  evidence  that  the  rate  of  the  induction  phase  reactions  of 
bulk  TNT  below  400 *C  is  dominated  by  the  formation  of  2,4-DNA  and 
compounds  closely  related  to  it.  Moreover,  there  is  virtually 
undisputable  augmenting  mechanistic  evidence  discussed  above  that 
both  intra-  and  intermolecular  H«  transfer  from  -CH3  to  -NOz 
initiate  thermal  decomposition  of  TNT  below  400 *C.  To  represent 
the  rate  of  this  induction  phase  process  in  bulk  TNT,  the  E.  and  A 
values  of  Cook  and  Abegg  (Reference  79)  are  useful  because  they 
are  about  the  average  of  the  data  available  (Table  5) . 

As  discussed  in  Section  IV. 1.,  C-N02  homolysis  is  the  dominant 
initial  reaction  of  nitrobenzene  compounds  above  about  900 ’C.  To 
assess  this  reaction  in  TNT,  the  rate  data  for  1,3,5-TNB  in  Table 
6  are  valuable.  Finally,  for  the  acceleratory  phase,  the  global 
rate  data  reflect  the  effects  of  catalysis  on  the  numerous 
reactions  that  take  place  in  this  stage.  The  dominant  reactions 
are  unknown,  but  catalysis  is  an  important  factor  in  the  rate. 

The  rate  data  of  Robertson  (Reference  44)  in  Table  5  are  useful 
for  representing  this  phase.  These  three  reactions  lead  to 
conclusion  (b) . 


(b)  The  proper  description  of  the  dominant  initiation  reaction 
of  bulk  TNT  depends  on  the  temperature  and  reaction  time.  This 
statement  follows  from  consideration  of  the  relative  importance  of 
the  three  classes  of  reactions  highlighted  in  (a)  in  terms  of  the 
time  required  to  complete  the  reaction  at  a  given  temperature. 

This  can  be  assessed  by  equation  (29) .  Figure  4  is  a  plot  of  the 
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(29) 


T  -  _ -2 _ 

R  (In A  ♦  InAt) 

results.  Note  that  the  rate  of  the  bulk  phase  Induction  reaction 
corresponding  to  the  formation  of  2,4-DNAn  and  the  oxidation  of 
-CH3  matches  the  rate  of  C-N02  homolysis  at  about  770*C.  Above 
770 *C,  C-N02  homolysis  is  the  faster  reaction,  while  below  770 *C, 
the  -CH3  oxidation  reaction  is  faster.  The  reaction  time  at  this 
crossover  temperature  is  about  l  msec.  Consequently,  the  impact 
sensitivity  which  results  in  heating  for  about  250  Msec  might  be 
expected  to  be  influenced  by  reactions  of  the  -CH3  group,  as 
Kamlet  and  Adolph  (Reference  24)  and  Bliss,  fit  al.  (Reference  181) 
found.  However,  for  shock  sensitivity,  in  which  the  full  energy 
input  occurs  in  about  1  Msec  or  less  and  T  >  1000 *C,  ono2 
homolysis  contributes  more  than  the  -CH3  group  oxidation.  It  is, 
therefore,  understandable  that  the  deuterium  isotope  effect 
(Reference  66)  is  diluted,  because  C-N02  homolysis  has  become  very 
important.  The  fact  that  the  impact  sensitivity,  Hso,  is  more 
influenced  by  -CH3  group  oxidation,  while  the  shock  sensitivity, 
P90,  is  more  influenced  by  C-N02  homolysis,  is  the  probable  reason 
why  TNT  is  an  outlying  point  on  a  plot  of  Hso  vs.  P90  for  a  series 
of  explosives  (Reference  27) . 

(c)  Times-to-explosion  of  TNT  longer  than  0,1  sec  are 
controlled  bv  catalytic  reactions  of  the  accelerators  phase. 
Shorter  .thermal  esplpgiQD  times  are  . were  controlled  by  c-nq2 
homolvsis  and  -CH-.  oxidation.  Experimental  times-to-explosion 
measured  in  different  ways  over  eight  decades  of  time  are  given 
for  TNT  in  Figure  4.  Notice  that  between  1  and  60000  sec  the 
time-to-explosion  is  not  modeled  by  the  induction  phase  kinetics 
of  Table  1,  as  might  be  expected,  but  rather  by  the  acceleratory 
phase  kinetics  which  reflects  catalysis.  At  shorter  reaction 
times  (10"2  -  10~3  sec)  the  times-to-explosion  depart  from  the 
acceleratory  phase  kinetics.  This  could  simply  be  the  result  of 
experimental  inadequacies,  such  as  the  inability  to  achieve 
efficient  heat  transfer,  that  place  the  measured  temperature 
higher  than  the  true  sample  temperature.  Acknowledging  the 
experimental  shortcomings,  the  times-to-explosion  in  Figure  4  at 
short  times  depart  toward  the  induction  puase  rate  and  the  C-N02 
homolysis  rate  and  away  from  the  catalytically  dominated 
acceleratory  phase  rate.  Hence,  at  short  reaction  times  primary 
reaction  steps  increasingly  determine  the  time-to-explosion. 
However,  at  longer  reaction  times,  catalysis  and  secondary 
reactions  dominate  the  time-to-explosion.  It  is  no  surprise, 
therefore,  that  E,  values  from  times-to-explosion  measurements 
have  never  been  correlatable  with  the  molecular  structure  and 
elementary  reaction  kinetics  of  polynitroaromatic  compounds. 
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Eq.  29  (C-N02  HOMOLYSIS 
Eq.  29  (INDUCTION) 

Eq.  29  (ACCELERATORY) 
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Figure  4.  Ti*es-to-Explosion  for  TNT  Modelled  with  Rate  Constants  of  Three 
Processes:  C-N02  Homolysis,  -CH,  Oxidation  (Induction  Phase),  and 
Catalysis  by  Decomposition  Products  (Acceleratory  Phase) . 


b.  o-HH  Bonds 


Ths  difficulty  of  correlating  ths  sxplosivs  characteristics 
of  fi-aainonitrobenzene  compounds  with  their  molecular  and  crystal 
features  was  described  in  Section  XX. 4.  The  amino  group 
stabilizes  the  nitroaromatic  molecule  toward  shock  sensitivity 
(Reference  67)  and  impact  sensitivity  (Reference  24,181),  but 
hexanitrobiphenyl  explosives  are  exceptions  in  that  the  impact 
sensitivity  is  actually  increased  by  the  presence  of  -NH2 
(Reference  183) .  Xn  these  latter  compounds,  steric  demands  may 
weaken  several  of  the  C-N02  bonds.  Xn  reality,  many  variables 
change  when  “11112  groups  are  added  to  the  trinitrobenzene  backbone 
making  it  unlikely  that  a  single  parameter  holds  the  key  to  the 
dramatic  changes  in  the  chemical  and  physical  properties  through 
the  series  TNB,  MATE,  DATB,  TATB.  The  differences  are  partly  due 
to  the  crystal  structure  and  partly  due  to  the  molecular  and 
electronic  structure.  The  important  features  of  the  crystal 
structure  are  considered  first. 

Structurally,  TNB,  MATB,  DATB  and  TATB  are  very  different. 

The  intermolecular  bond  distances  were  given  in  Table  3.  TNB 
contains  both  planar  and  non-planar  benzene  rings  with  -N02  groups 
twisted  and  bent  out  the  plane  and  0  atoms  with  large  thermal 
parameters  (Reference  68) .  The  -N02  groups  of  MATB  are  non- 
coplanar  with  the  benzene  ring  (Reference  69) .  The  dihedral 
angles  are  22.5*,  4.0*  and  8.5*.  Curiously,  the  most  rotated  -N02 
group  is  2-N02  which  might  have  been  expected  to  be  coplanar 
because  of  intramolecular  hydrogen  bonding.  Instead, 
intermolecular  hydrogen  bonds  force  this  -NOz  group  out  of  the 
plane.  From  the  point  of  view  of  the  structure  alone  DATB 
(Reference  70)  and  TATB  (Reference  71)  are  relatively  similar. 

Both  molecules  are  planar  with  chain-like  intermolecular  hydrogen 
bonding  in  DATB  and  layer-like  bonding  in  TATB. 

Electronically,  the  presence  of  the  -NHZ  and  -N02 
substituents  on  an  aromatic  ring  produces  intricate  competitive 
effects.  The  -NH2  group  is  inductively  withdrawing  but  resonance 
donating  of  electron  density,  whereas  the  -NOz  group  is  both 
resonance  and  inductively  withdrawing.  Therefore,  it  is  difficult 
to  predict  the  reactivity  of  the  ring  when  both  groups  are 
present.  Politzer,  al.  (Reference  184)  and  Rogers,  al. 
(Reference  67)  addressed  this  problem  computationally.  The 
addition  of  an  -NH2  group  in  the  ortho  position  of  nitrobenzene 
stabilizes  the  ring  because  of  resonance.  Accordingly,  the  C-N02 
bonds  become  progressively  shorter  as  the  number  of  -NH2  groups 
increases  (Table  3) .  SCF-3-21G  MO  calculations  (Reference  185) 

reveal  no  energy  minimum  for  the  aci  tautomer  (VI)  compared  to  o- 


(VI) 


nitroaniline.  However,  a  significant  body  of  experimental  data  is 
at  odds  with  these  computational  findings.  First,  Matveev,  fi£  al. 
(Reference  102)  observed  a  primary  deuterium  isotope  effect  (Table 
8)  in  the  gas  phase  when  o-nitroaniline  and  o-nitroaniline-d2 
thermally  decompose.  This  indicates  that  N-H  bond  scission, 
possibly  of  the  type  shown  in  (VI),  occurs  in  the  RDS.  E.  for  the 
decomposition  step  is  relatively  large  (57.5  kcal/mol)  which 
necessitates  a  relatively  high  temperature  for  the  reaction  to 
occur  at  a  rapid  rate.  This  temperature  is  still  below  that  of 
other  plausible  reaction  routes,  such  as  C-N02  homolysis  for  which 
the  E.  values  are  65  -  70  kcal/mol  (Tables  6  and  7) .  Second,  the 
principal  mass  spectral  fragmentation  step  of  ortho- 
aminonitrobenzene  compounds  is  loss  of  OH  (References  36,126). 

This  fact  strongly  suggests  that  H  transfer  as  in  (VI)  is  facile. 
Third,  ESR  spectra  of  irradiated  TATB  produces  long-lived  radicals 
consistent  with  protonated  TATB  (VII) ,  which  is  proposed  to  result 


(VII) 


from  intermolecular  H*  transfer  (References  186,187).  In  view  of 
positive  evidence  for  species  or  processes  involving  or  formally 
resembling  (VI)  in  either  the  transition  state  or  as  the  stable 
form,  the  negative  computational  evidence  (Reference  185) ,  evokes 
skepticism.  More  consistent  with  experiment  is  the  computational 
finding  that  conversion  of  2-nitroaniline  to  benzofurazan  by 
equation  (14a)  is  nearly  thermally  neutral  (Reference  152) .  There 
is  substantial  experimental  evidence  discussed  below  for  these 
conversions  when  -NH2  and  -N0Z  groups  are  ortho  related. 

Equation  (14a)  has  been  known  for  at  least  80  years  and  is 
especially  prevalent  in  an  oxidizing  environment  (References 
188,189).  These  reactions  are  standards  in  organic  synthesis 
(References  141,190).  In  addition  to  the  primary  deuterium 
isotope  effect  when  2-nitroaniline  decomposes  (Reference  102) , 
Rogers,  fit  al .  (Reference  191)  found  a  primarily  isotope  effect 
(kg/ko  *  1.5)  when  TATB  thermally  decomposes  at  297  -  382 *C. 

These  findings  support  equation  (14)  as  the  lower  temperature 
decomposition  channel.  The  unimolecular  decomposition  of  MATB  is 
a  straightforward  extension  of  the  behavior  of  2-nitroaniline  in 
that  it  loses  OH  by  MS/MS-CID  (Reference  126) .  However,  in 
addition  to  the  potential  for  cyclization  reactions,  equation 
(14) ,  MATB  also  loses  NO  from  the  radical  cation  produced  by  CID 
(Reference  126) .  The  mechanism  is  proposed  to  be  equation  (30) , 
whose  product  is  resonance  stabilized.  Furazan  and  furoxan 
formation  occurs  with  TATB  that  has  been  thermally  decomposed  and 


subcritically  shocked  or  impacted  (s*f«rmc«a  162,192-194).  At 
least  thy  different  furaiani  were  wr.x  ctid  by  TLC  and  analysed 


by  CI-MS  whan  TATB  was  tharaally  dacoaposad  (Rafaranca  192) . 

Curva  raaolution  of  tha  XPS  1 Inashapaa  of  TATB  aub j  actad  to  shock 
and  ispact  ara  conalstant  with  formation  of  about  6%  furazan  and 
6%  furoxan  darivativas  (Rafaranca  193) .  On  tha  othar  hand,  there 
is  no  ZR  evidence  of  tha  furazan  ring  in  thermally  decomposed  TATB 
(Rafaranca  195) . 


Kinetics  information  on  MATB,  DATB,  and  TATB  in  the  solid, 
liquid  and  gas  phase  is  compiled  in  Table  10.  Data  for  several 
alkylamine  derivatives  of  2,4,6-trinitrobenzene  are  not  shown  but 
are  available  (Reference  54) .  The  tendency  for  these  compounds  to 
react  by  aquation  (14)  is  tha  most  logical  explanation  for  the 
much  smaller  values  of  Ea  and  A  compared  to  TNB  (Table  6) . 
Unfortunately,  further  analysis  is  necessarily  qualitative  because 
of  scatter  in  the  data  and  the  tact  that  the  measurements  were 
made  by  different  methods  (DTA,  DSC,  manometry,  electron  beam 
heating  traces,  and  affusion  cell  mass  spectrometry) .  In  the 
liquid  and  solid  phase,  matrix  effects  and  catalysis  by  the 
decomposition  products  influence  the  Arrhenius  constants  (Section 
II. 3.).  There  is  a  rough  pattern  of  increasing  E,  and  A  from  MATB 
to  TATB,  which  is  consistent  with  the  decrease  in  impact 
sensitivity  along  this  series  (Table  1) .  However,  the  Arrhenius 
factors  of  TNB  are  clearly  out  of  line  because  the  initiation 
process  (C-M02  homolysis)  is  different  from  equation  (14) .  To 
analyze  the  relative  roles  of  C-N02  homolysis  and  furazan/ furoxan 
cyclization  in  the  aminonitrobenzene  series,  the  gas  phase  data 
for  MATB  in  Table  10  are  perhaps  least  encumbered  by  matrix 
effects  and  catalysis.  These  rate  data  are  compared  in  Figure  3 
to  that  of  TNB  in  Table  6  (Reference  198) .  The  temperature  at 
which  C-N02  hemolysis  equals  cyclization  is  475 *C  or  625 *C 
depending  on  the  MATB  rate  chosen.  Of  course,  these  are  crossover 
temperatures  for  the  gas  phase.  Their  relevance  to  the  condensed 
phase  is  uncertain.  If  trusted,  these  temperatures  are  easily 
achieved  in  the  condensed  phase  during  shock  and  impact  where 
Wenograd  (Reference  2)  estimated  that  TNB  reaches  1060 *C  in  about 
250  Msec.  Estimates  of  the  crossover  temperatures  using  the  other 
data  in  Table  10  for  MATB,  DATB  and  TATB  produce  unreasonably 
large  temperatures  or  impossible  negative  temperatures . 
Consequently,  it  is  uncertain  at  which  temperature  C-N02  homolysis 
becomes  competitive  with  cyclization  in  these  compounds.  However, 
Owens,  fit  al.  (References  73,74)  and  Murray,  At  si.  (Reference  75) 
proposed  that  characteristics  of  the  C-N02  bond  are  the  main 


Table  10 


Arrhaniua  Data  for  Thermal  Decomposition  of  Amino- 
2,4,  6-Trinitrobenzene  Compounds 


Compound 

Phase 

T,  *C 

Ea,  kcal/mol 

A,  sec*1 

Ref. 

MATB 

gas 

301  -  343 

38.5 

10#i 

53 

gas 

240  -  280 

41.8 

1011 

84 

liquid 

240  -  280 

32 

10,a 

84 

liquid 

249  -  259 

48 

- 

54 

DATB 

liquid 

320  -  323 

46.3 

1013.07 

12 

liquid 

- 

50 

1015 

58 

solid 

220  -  270 

47 

1013.2 

58 

TATB 

gas/solid 

250  -  300 

43 

- 

196 

solid 

297  -  382 

59.1 

101«37 

191 

solid 

331  -  332 

59.9 

1Q195 

12 

solid 

345  -  375 

62±4 

1019 

197 

51 


or  rate  controlling  factor  in  the  shock  and  impact  sensitivity. 
This  is  the  only  initiation  reaction  that  TNB,  MATB,  DATB,  and 
TATB  could  have  in  common  to  account  for  the  smoothly  increasing 
trend  of  impact  sensitivity  (Table  1) .  The  observation  of  furazan 
and  furoxan  compounds  in  heated,  impacted  and  shocked  samples  of 
TATB  (References  162,192-194)  is  understandable  from  the  fact  that 
the  temperature  during  heat-up  and  cool-down  emphasizes  the  lower 
Ea  cyclization  channel.  The  shocked  and  impacted  samples  did  not 
explode  and,  therefore,  did  not  reach  the  temperature  range  where 
C-N02  homolysis  could  dominate.  Essentially  the  same  spectral 
features  subsequently  interpreted  as  furazan  and  furoxan  formation 
(References  162,192-194)  were  once  interpreted  as  evidence  of  C- 
N02  homolysis  (References  119,200). 

C-N02  homolysis  lacks  appeal  as  the  source  of  the  decrease  in 
shock  and  impact  sensitivity  as  -NH2  groups  are  added.  This  is 
because  the  C-N02  BDE  is  not  especially  sensitive  to  the  other 
substituents  on  the  ring  (References  32,52).  Thus,  while  C-N02 
homolysis  is  the  only  common  chemical  initiation  process  possible 
TNB,  MATB,  OATB,  and  TATB,  the  increase  in  the  amount  of  inter- 
and  intramolecular  hydrogen  bonding  along  this  series  is  a 
structural  variable  that  may  influence  the  shock  sensitivity. 

These  interactions  have  been  proposed  by  Rogers,  gt  al.  (Reference 
67)  to  absorb  shock  energy  and  prevent  it  from  being  localized  in 
the  molecule.  This  effect  is  difficult  to  judge  kinetically,  but 
it  could  be  accommodated  in  theories  of  how  mechanical  energy 
converts  into  chemical  processes  (References  14-19) .  Raman 
spectroscopy  during  dynamic  shock  loading  also  suggests  that  the 
N-H* • • O  bond  becomes  compressed  (Reference  35).  A  qualitative 
indication  of  the  increase  in  the  cohesive  force  in  the  lattice  as 
-NH2  groups  are  added  is  the  increase  in  melting  points  (Table  2) . 

Taken  together,  the  cyclization  reactions  (equation  14)  of 
aminonitrobenzene  explosives  dominate  the  initial  thermal 
reactions  in  the  lower  temperature  range.  These  reactions  could 
be  important  in  sensitizing  a  "rough  handled"  explosive  (Reference 
78) .  At  high  temperatures  characteristic  of  shock  and  impact 
initiation,  C-N02  homolysis  appears  to  play  a  large  and  perhaps 
even  dominant  role  in  the  initiation  chemistry.  However,  the 
decreased  shock  and  impact  sensitivity  as  -NH2  groups  are  added  to 
the  ring  seems  to  be  best  accounted  for  by  shock-absorbing 
hydrogen  bonding  in  the  crystal  lattice. 

TATB  is  the  most  important  explosive  of  this  category  and,  as 
such,  deserves  further  comment.  TATB  is  infamous  for  its  high 
thermal  stability  (>300*C),  low  impact  sensitivity,  and  low 
solubility.  The  C-NH2  bonds  are  unusually  short  (1.319(7)  A), 
while  the  C-C  bonds  are  unusually  long  (1.444(7)  A)  compared  to  C- 
C  of  benzene  (1.40  A)  (Reference  71).  A  typical  C-C  single  bond 
distance  is  about  1.47  A.  The  n*  MO's  of  the  ring  are  partially 
occupied  (Reference  67).  Strong  intramolecular  NH* • • O  hydrogen 
bonding  leads  to  a  planar  molecule  while  moderate  intermolecular 
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hydrogen  bonding  causes  stacking  of  the  molecules  in  graphite-like 
layers  (References  71,201). 

Given  the  unusual  materials  properties  of  TATB  and  the 
significant  amount  of  intermolecular  association,  complexity  in 
the  decomposition  behavior  of  TATB  comes  as  no  surprise.  Unlike 
most  nitroaromatic  explosives,  TATB  decomposes  from  the  solid 
phase  unless  it  is  heated  rapidly  (Reference  197).  A  melting 
point  of  about  450*0  has  been  mentioned  (Reference  197)  although 
Meyer  (Reference  13)  gives  350* C.  Discoloration  rather  than 
melting  is  visually  observed  when  TATB  is  heated  at  2000*0/ sec  to 
500*0  on  a  Pt  filament  (Reference  22).  Catalano,  al. 
(References  45,195,202)  made  numerous  observations  and 
measurements  on  the  slow  thermal  decomposition  of  TATB.  Reactions 
resulting  in  a  residue  and  gas  formation  occur  simultaneously  with 
sublimation  in  the  305  -  357*0  range  (Reference  45).  By  DSC,  the 
process  is  exothermic  despite  the  fact  that  AH^p  is  endothermic  by 
40  keal/mol  (References  203,204).  Perhaps  five  processes  occur  at 
342  -  357*0  necessitating  the  use  of  multiple  rate  expressions 
(Reference  45) .  The  initial  stage  is  autocatalytic  with  E.  ■  50  - 
60  keal/mol  and  appears  to  consist  of  at  least  two  consecutive 
reactions.  The  global  reaction  order  is  0.5.  The  overall  Ea 
resembles  the  values  given  in  Table  10.  The  remaining  stages  of 
decomposition  have  Ea  values  in  the  range  of  50  -  80  keal/mol 
(Reference  3) . 

MS  and  IR  analyses  of  the  gaseous  products  from  confined  and 
unconfined  TATB  suggest  (Reference  202)  that  the  removal  of 
intercalated  C02,  acetone,  H20  and  N2  occurs  initially.  1  -  2%  of 
the  space  between  the  layers  is  occupied  by  these  trapped 
molecules.  It  is  conceivable  that  the  small  amount  of  H2  also 
observed  comes  from  furazan  formation,  equation  (14a),  while  some 
of  the  H20  may  result  from  furoxan  formation,  equation  (14b) .  As 
the  reaction  rate  accelerates  into  an  explosion  C02,  H20  and  N2 
increase  in  concentration  and  are  joined  by  HCN,  CO  and  NO.  Other 
gaseous  products  detected  in  small  amounts  included  N20,  N02,  CH* 
and  C2N2.  A  solid  residue  containing  C,  H,  N,  and  0  remained. 

The  gas  products  from  rapid  thermolysis  of  TATB  are  available  in 
other  studies  (Reference  198)  and  are  discussed  in  Section  VI. 
Furazans  and  furoxans  are  detected  in  some  thermal  decomposition 
studies  of  TATB  (References  162,192-194)  but  not  others, 

(Reference  195) .  These  compounds  are  not  likely  to  have  a  high 
steady-state  concentration  because  they  are  less  thermally  stable 
than  TATB.  The  gas  products  of  decomposition  of  furoxans  closely 
resemble  those  of  TATB  (Reference  205) .  It  should  be  noted  that  a 
completely  different  view  of  thermal  decomposition  of  TATB  has 
been  presented  from  MS  studies  by  Farber  and  Srivastava  (Reference 
196) .  The  El  mass  spectrum  of  the  products  following  use  of  an 
effusion  cell  and  Langmuir  evaporation  at  200  -  300 °C  showed 
sublimation  followed  by  splitting  of  TATB  into  two  fragments  of 
m/e  *  114  and  144.  These  two  masses  are  consistent  with  fission 
of  the  molecule  along  the  dashed  line  of  (VIII) .  Species  having 
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other  m/m  values  mxm  detected  as  wall  that  would  rssult  fro* 
various  sntitias  claavad  diraetly  fro*  TATB.  Thasa  rasults  ara 
uniqua  and  parhaps  raflact  tha  conditions  of  tha  experiment  aora 


(VIII) 


than  tha  relevance  to  tha  bahavior  of  TATB  as  a  bulk  axplosiva. 

For  instanca ,  thay  contradict  tha  primary  deuterium  isotope  effect 
found  during  decomposition  of  bulk  TATB  (Reference  191) .  The 
value  of  E.  -  43  kcal/mol  in  Table  10  from  MS  studies  of  TATB 
(Reference  196)  was  determined  from  the  rate  of  growth  of  m/e  * 
114.  This  value  of  E*  is  disquietingly  similar  to  AH^  of  TATB  of 
40  kcal/mol.  Moreover ,  Langmuir  evaporation  and  effusion  cell 
rasults  for  another  molecule,  HMX,  sharply  differ  from  those  of 
another  MS  technique ,  simultaneous  thermogravimetry-modulated  beam 
mass  spectrometry  (Reference  38) . 

c.  a -OH  Bonds 

The  -OH  group  plays  a  major  role  in  the  thermal  stability  of 
nitroaromatic  compounds ,  but  the  effect  has  to  do  with  the 
tendency  of  -OH  to  become  a  quinoid  group,  rather  than  with 
cyclization.  The  -OH  group  is  a  resonance  electron  donor  as  is 
the  -NHa  group.  However,  the  -OH  substituent  shortens  the  C-N02 
bond  and  lengthens  the  N-0  bond  less  than  does  -NH2  (References 
185,206).  Since  -NH2  stabilizes  the  nitrobenzene  backbone  toward 
initiation  by  impact  (Table  1) ,  it  might  be  expected  that  -OH  also 
stabilizes  the  molecule  compared  to  TNB.  The  opposite  is  found  in 
practice  because  -OH  derivatives  are  more  impact  sensitive  than 
TNB  (Reference  24) .  This  has  recently  been  explained  by  the  fact 
that  the  aci-nitro  structure  (X)  is  computed  to  have  an  energy 
minimum,  but  at  only  1.54  kcal/mol  above  (IX)  (Reference  185). 

Upon  shock  or  impact  stimulus,  vibrational  excitation  could  force 


a  significant  fraction  of  the  molecule  into  the  more  reactive  form 
(X) .  This  computational  finding  (Reference  185)  is  consistent 
with  a  century  of  observations  in  organic  chemistry.  Quinones  are 
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isolated  upon  thermal  decomposition  of  p-nitrophenols ,  equation 
(31)  (References  207,208).  The  quinoid  fora  (X)  has  long  been 
suspected  to  be  the  source  of  bright  color  when  picric  acid  (PA) 
is  dissolved  in  basic  solution.  For  this  reason,  nitrophenol 
coapounds  are  useful  as  acid-base  indicators. 


PA  is  a  historically  important  explosive,  but  is  hazardous 
because  of  its  high  sensitivity  to  shock.  The  acidity  of  PA  (k.  - 
1.6  x  10'1  M)  is  also  detracting  because  it  corrodes  aetal  and 
attacks  aany  organic  aaterials  with  which  PA  might  be  mixed. 
Despite  these  problems,  PA  was  used  extensively  as  an  explosive  up 
to  World  War  II. 

Compared  to  the  -CH3  and  -NH2  groups  relatively  little  work 
has  been  conducted  on  the  details  of  thermal  decomposition  of 
hydroxy-substituted  nitrobenzene  compounds  (Reference  209) .  The 
behavior  of  PA  is  reminiscent  of  TNT  and  other  nitroaromatics  in 
that  melting  followed  by  decoaposition  occurs.  The  process  is 
first-order  &utocatalytic  after  the  initial  induction  period 
(Reference  209) .  Like  TNT,  the  rate  of  decomposition  in  the  gas 
phase  is  slower  than  in  the  liquid  state  suggesting  that  catalysis 
and  intermolecular  interactions  are  important.  Arrhenius  data  are 
sparse,  but  are  shown  in  Table  11.  The  large  variation  in  values 
makes  conclusions  difficult  to  draw.  However,  the  TTX  value 
(Reference  94)  probably  can  be  discarded  because  of  the  difficulty 
of  obtaining  efficient  heat  transfer. 

The  initial  reaction  step  of  PA  at  the  temperatures  in  Table 
11  appears  to  involve  the  quinoid  structure  (X)  on  the  basis  of 
ESR  and  MS  studies  (References  126,132,170,211).  Radicals  (XI) 
and  (XII)  that  are  readily  formed  from  the  equivalent  of  (X)  were 
suggested  to  explain  the  ESR  spectra  of  thermally  decomposing  PA 
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Table  11.  Arrhenius  Data  for  Hydroxy- substituted 
Trinitrobenzene  Compounds 


Compound 

T,  *C 

Et,  keal/mol 

A,  sec*1 

Method 

Ref 

PA 

183  -  270 

38.6 

1QU.7 

EG 

210 

233  -  244 

51.9 

- 

DTA 

56 

267  -  350 

27.4 

- 

TTX 

94 

1  #  3- (OH)  2-2- 
4 , 6-TNB 

180  -  200 

34.6 

IqU.2 

EG 

2 10 

1,3,5-<0H)S- 
2,4, 6-TNB 

145  -  165 

54.3 

1022 

EG 

50 

56 


(Reference  211) .  A  thermal  decomposition  reaction  that  appears  to 
be  unique  to  PA  (and  MATB)  is  inferred  from  the  ions  in  Cl -MS 
(Reference  170)  and  MS/MS-CID  (References  126,132).  This  is  loss 
of  NO  as  the  dominant  decomposition  pathway.  In  a  manner  similar 
to  equation  (31) ,  resonance  stabilization  through  the  quinoid 
structure  according  to  equation  (32)  was  proposed  to  be  the  main 
reason  for  the  facileness  of  this  process.  When  all  data  are 
taken  together,  a  clear  understanding  of  the  decomposition 


mechanism  of  PA  does  not  exist.-  The  mobility  of  the  H  atom  bound 
to  -OH  almost  certainly  plays  a  major  role  in  thermal 
decomposition  mechanism  at  lower  temperatures.  Unlike  when  -CH3 
and  -NH2  groups  are  ortho  to  -N02,  cyclization  does  not  occur.  It 
is  logical  to  extrapolate  from  the  more  extensive  data  available 
on  — CH3  and  -NH2  substituted  trinitrobenzene  compounds  that  C-N02 
homolysis  is  also  a  major,  if  not  the  dominating,  early  reaction 
during  impact  and  shock  initiation  of  PA. 


d.  -n3  ortho  to  -no2 


The  azide  group  cyclizes  an  ortho  -N02  group  in  a  manner 
resembling  the  behavior  when  -CH3  and  -NH2  groups  are  ortho  to 
-N02.  Upon  heating  above  70 *C  2-nitrophenylazide  eliminates  N2 
and  forms  benzofuroxan  according  to  equation  (33).  This  reaction 
has  been  used  for  years  as  a  facile  synthetic  route  to  substituted 
benzofuroxans  (References  141,190,212).  In  the  El  mass  spectrum 


of  2-nitrophenylazide,  an  m/e  =  136  fragment  is  observed 
corresponding  to  benzofuroxan  (Reference  64).  Equation  (33) 
contrasts  with  the  behavior  of  meta  and  para  isomers  where  the 
rate  of  thermolysis  is  10*  times  slower  at  100 *C  owing  to  a 
different  decomposition  mechanism,  equation  (34)  (References 
106,213) . 
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The  Arrhenius  data  compiled  in  Table  12  reflect  the 
significantly  different  rates  of  equations  (33)  and  (34).  When 
ortho  related  -N3  and  -N02  groups  are  present  on  the  ring,  the 
decomposition  rate  is  insensitive  to  the  presence  of  other  -N3 
or  -H02  groups,  because  equation  (33)  dominates  the  process 
(Reference  106).  On  the  other  hand,  equation  (33)  is  accelerated 
by  the  presence  of  electron-withdrawing  groups  p&ca  to  -N02 
(Reference  215) . 

Whether  equation  (33)  or  (34)  dominates  has  a  pronounced 
effect  on  the  burn  rate  of  these  compounds.  Despite  the  higher  Ea 
of  equation  (34) ,  the  bum  rates  of  3-  and  4-nitrophenylazide  are 
1.4  -  1.7  times  higher  than  2-nitrophenylazide  (Reference  106). 
This  is  because  the  nitrene  product  of  equation  (34)  decomposes 
much  more  exothermically  than  the  furoxan  product  of  equation 
(33) .  Release  of  a  large  amount  of  heat  early  in  the 
decomposition  reaction  scheme  near  the  burning  surface  produces  a 
high  bum  rate  for  the  material  (References  216,217). 

The  initial  decomposition  chemistry  when  an  azide  group  is 
present  on  a  nitroaromatic  backbone  appears  to  be  straight¬ 
forward  based  on  current  knowledge  of  these  compounds.  However, 
the  explosive  characteristics  of  these  compounds  are  essentially 
unknown. 

e.  Other  ortho  Substituents 

Cyclization  of  an  -N02  group  with  an  ortho  substituent  is  so 
widely  occurring  that  a  few  additional  studies  are  worthy  of 
mention  to  dispel  any  lingering  thought  that  the  results  in 
Sections  IV. 4. a,  b,  and  d  are  unique. 

As  noted  in  Section  IV. 4. a,  a-CH  abstraction  from  a  -CR2 ' R" 
group  (R'  *  R"  *  H?  R'  *  H,  R"  -  CH3;  R‘  =  CH3,  R"  =  H)  is  a 
general  reaction  in  the  presence  of  an  ortho  -N02  group.  However, 
further  reaction  of  the  molecule  following  this  initial  C-H 
abstraction  step  proceeds  by  different  routes  (Reference  37) .  For 
instance,  the  formation  of  anthranil  by  equation  (26)  occurs  for 
the  -CH3  group.  On  the  other  hand,  2-nitroethylbenzene  yields 


(XHD 


indole  (XIII)  through  a  series  of  H  rearrangements  (Reference  37) 
When  R'  *  H  and  R"  =  N(CH3)2  equation  (35)  occurs  in  the 
sample  on  sitting  at  room  temperature  for  several  days  (Reference 
218) . 
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Table  12.  Arrhenius  Data  for  Thermal  Decomposition  of 
Nitrophenylazide  Compounds 


Compound 

Ea,  Kca  1/mol 

A,  sec*1 

Ref 

2 -Nitrophenylazide 

26.1 

1012.3 

214 

4 -Nitrophenylaz ide 

40.6 

10168 

213 

1,3, 5-Tr iaz ido-2 ,4,6 

26.0 

1012.1 

106 

-trinitrobenzene 
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CX* 


-  H ,  "  CHj 

-20 


Unimolecular  decomposition  of  2 -nitroanisole  has  bean  proposed 
to  occur  via  M  migration  and  ring  closure,  equation  (36) 

(Reference  144) .  This  species  then  collapses  to  benzoxazole 
(Reference  37) . 


Such  a  rearrangement  is  consistent  with  the  formation  of  a 
substituted  benzoxazole  when  a  longer  alkoxy  chain  is  present, 
equation  (37)  (Reference  219).  Likewise,  2-nitro(methylamino) 


benzene  was  observed  to  form  benzimidazole  at  a  relatively 
relatively  low  temperature  by  equation  (38)  (Reference  37) . 
Benzofurazan  is  isolated  from  methyl-N-(2-nitroaryl) carbamate 

(38) 

according  to  equation  (39)  (References  220,221). 


or" 


-2H 


H 


N0a 


-2  0 


CO* 


^^^NHCOaCH*  250-270°C 


The  hydrazine  derivative  of  nitrobenzene  also  decomposes  at  a 
relatively  low  temperature  liberating  water  and  forming 
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bansotriazola,  aquation  (40)  (Rafaranca  123).  Evan  1,2- 
dinitrobans ana  losas  0a  upon  haatinq  giving  snail  amounts 


of  banzofurazan  and  banzofuroxan  (Rafarancas  102,110)  by  aquation 
(41) . 
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SECTION  V 


EFFECT  OF  SUBSTITUENTS  ON  THE  KINETICS  OF  THERMOLYSIS 

Throughout  Section  IV,  the  positioning  of  a  reactive 
substituent  ortho  to  the  -N02  group  stands  out  as  a  very  important 
factor  in  the  low  Ea  thermal  decomposition  pathways  of 
nitroaromatic  compounds.  In  this  section  a  broader  view  is  taken 
of  the  effect  of  various  substituents  on  the  decomposition  rate. 

Hammett  championed  the  notion  that  rate  constants  of  various 
reactions  are  related  to  the  structure  and  composition  of  the 
parent  aromatic  molecule  (References  222,223).  More  specifically, 
the  electron  density  of  the  aromatic  ring  influences  its 
reactivity  and  the  reactivity  of  the  attached  substituents. 

Linear  free  energy  relationships  for  the  substituents  were 
determined  that  are  frequently  referred  to  as  Hammett  or  Taft 
constants  (References  224-226) . 

An  outstandingly  prominent  feature  of  polynitroaromatic 
compounds  is  that  thermal  decomposition  data  generally  do  not 
correlate  with  Taft  or  Hammett  constants.  The  absence  of 
correlation  is  independent  of  the  phase  of  decomposition  in  that 
neither  the  gas  phase  nor  the  condensed  phase  kinetic  data 
correlate  (Reference  52).  Several  possible  reasons  for  non¬ 
correlation  are  possible.  Hammett  constants  depend  on  the  solvent 
and  the  temperature  of  measurement.  First,  the  temperature  at 
which  polynitroaromatic  compounds  decompose  is  usually  very 
different  from  the  temperature  at  which  the  Hammett  constants  are 
determined.  Second,  the  molecular  environment  during  thermal 
decomposition  is  different  from  that  in  which  the  Hammett 
constants  are  determined.  Third,  the  foregoing  sections 
illustrate  that  catalysis  and  bimolecular  processes  frequently 
dominate  during  thermal  decomposition  of  nitroaromatic  compounds. 
Fourth,  the  rates  of  thermal  decomposition  of  nitroaromatic 
compounds  may  not  depend  very  much  on  the  resonance  and  inductive 
electronic  effects  of  the  ,  >bstituents.  For  example,  E.  for 
thermal  decomposition  of  it robenzene  compared  to  a  variety  of 
1,3-  and  1,4-substituted  robenzene  compounds  differs  by  only  ±2 
kcal/mol  and  the  A  values  differ  by  only  a  factor  of  five 
(Reference  52) . 

Despite  the  absence  of  correlations  with  Hammett  constants, 
the  substituents  on  the  ring  still  exert  a  profound  effect  on  the 
explosive  properties  and  the  kinetics  of  decomposition.  The  most 
ambitious  attempt  to  discover  structure/reactivity/detonation 
relationships  has  been  the  work  of  Zeman,  s£  al*  (References  5-8, 
54).  That  a  relation  might  even  exist  could  be  traced  to  comments 
of  Cook,  si  al.  (References  3,4)  who  proposed  that  the  early 
decomposition  kinetics  affect  the  detonation  velocity  of  an 
explosive.  Zeman,  si  al. ,  based  their  conclusions  on  values  E. 
before  the  onset  of  autocatalysis;  TD,  the  temperature  immediately 
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after  tha  onset  of  the  exotherm;  P,  the  C-J  pressure;  D,  the 
detonation  velocity;  and  Q,  the  heat  of  explosion.  Equation  (42) 
was  used.  The  values  of  the  constants  a  and  b  depend  on  the 

log  (E/TD )  -  — L_5  (42) 

variable  P,  D2  or  Q  chosen.  By  using  data  for  74  compounds,  41  of 
which  are  nitroaromatic  compounds  without  other  energetic 
substituents,  Zeman  found  groupings  of  3  -  5  compounds  that  gave 
reasonably  linear  solutions  for  equation  (42).  One  group 
contained  10  compounds.  On  this  basis,  structure  and  electronic 
relationships  were  proposed  to  exist  between  the  early  kinetics  of 
slow  thermal  decomposition  and  detonation  (References  6,7).  Few 
of  these  groupings,  however,  are  chemically  satisfying.  For 
example,  the  largest  group  containing  10  compounds  contains  both 
sensitizing  (-OH)  and  desensitizing  (~NH2)  substituents  (Reference 
7) .  Another  contains  both  -CH3  and  -Cl  substituents  (Reference  5) 
which  do  not  decompose  by  the  same  mechanism.  In  fact,  Rogers,  et 
al.  (Reference  191)  would  dispute  the  foundation  of  Zeman 's 
assumed  connection,  at  least  for  TATB.  Zeman  (Reference  5)  also 
related  the  equilibrium  constant  for  the  activated  complex,  K*,  to 
D  on  a  log-log  plot  and  found  selected  groupings  of  compounds  in 
the  74  compound  data  base.  A  caveat  is  that  a  log-log  plot  tends 
to  be  relatively  insensitive  to  curvature.  Moreover,  some  of  the 
correlations  are  rather  loose,  which  partly  reflects  the 
difficulty  of  obtaining  accurate  experimental  data.  Zeman' s  DTA 
data  are  sensitive  to  the  sample  weight  (References  55,56,97-101). 

Considering  all  of  the  literature,  the  decomposition  and 
explosive  properties  of  nitroaromatic  compounds  are  clearly 
sensitive  to  the  substituents.  However,  most  relationships  defy 
simple  chemical  interpretation.  The  trends  in  the  liquid  phase 
and  the  gas  phase  are  usually  different  because  of  the  role  of 
catalysis,  bimolecular  processes,  and  the  environment. 
Nevertheless,  a  few  correlations  merit  further  mention. 

According  to  the  Hammett  and  Taft  constants,  the  -N02 
substituent  is  one  of  the  strongest  electron  withdrawing  groups 
available.  Accordingly,  the  ionization  potential  of  the  jt-HOMO  of 
the  benzene  ring  increases  as  the  number  of  -N02  groups  is 
increased  (Reference  227) .  Perhaps  for  this  reason,  the 
decomposition  rate  is  faster  as  the  number  of  -N02  groups  on  the 
ring  is  increased  (Reference  55) .  The  relative  rate  constants  at 
330*C  in  the  gas  phase  (Reference  53)  for  C6HsN02:  1, 3-(N02)2C6H4: 

1.3.5-  (N02)  3C6H3  are  1:20:34.  This  is  also  the  trend  in  the  bum 
rate  of  these  compounds  at  a  given  pressure  and  value  of  the 
pressure  exponent  (References  106,228).  The  position  of  the  -N0Z 
groups  relative  to  one  another  on  the  ring  affects  the  reaction 
rate  at  330*C.  This  is  evident  for  1,2,3-TNB,  1,2,4-TNB,  and 

1.3.5- TNB  which  have  relative  rates  of  1:0.175:0.03  (Reference 
53) .  This  trend  suggests  that  within  a  given  stoichiometry  the 
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thermal  stability  is  improved  by  making  the  molecule  more 
symmetrical  (Reference  229) .  Bliss,  al .  (Reference  181) 
explained  this  same  pattern  in  the  sensitivity  to  impact  in  terns 
of  the  greater  ease  of  nitro-nitrite  isomerization  when  -N0a 
groups  are  ortho  related. 

As  substituents  other  than  -N02  are  added  to  the  ring, 
resonance  and  inductive  effects  shift  the  f-MO's  relative  to  one 
another  in  complex  ways.  The  resulting  effect  on  the  rates  of 
thermal  decomposition  frequently  makes  random  chemical  sense.  For 
example,  -OH  is  thermally  destabilizing.  The  addition  of  -OH 
increases  the  decomposition  rate  in  the  order  (Reference  229) 

2,4, 6-trinitroresorsinol  >  2 , 4 , 6-trinitrophenol  >  1,3,5-TNB.  This 
additivity  effect  is  satisfying.  On  the  other  hand,  the  trend  for 
multiple  -CH3  groups  in  the  liquid  phase  is  not  additive.  One 
-CH3  group  destabilizes  the  molecule  while  two  and  three  -CH3 
groups  stabilize  it  relative  to  TNB;  e.o.  the  trend  in  thermal 
stability  is  TNT  <  TNB  <  2 , 4 , 6-trinitroxylene  <  2,4,6- 
trinitromesitylene  (Reference  229).  Other  evidence  exists  that 
the  effect  of  the  -CH3  group  on  a  nitroaromatic  backbone  is 
complicated.  For  example,  in  the  gas  phase,  the  2 -NT  is  less 
thermally  stable  than  nitrobenzene  because  of  anthranil  formation 
(Reference  33)  (Section  IV. 4. a),  while  1,4-NT  is  more  stable 
because  of  the  lower  A  factor  of  decomposition  (Reference  32). 

A  -Cl  substituent  destabilizes  2 , 4 , 6-trinitrochlorobenzene  as 
revealed  by  the  increase  of  1.75  in  the  gas  phase  reaction  rate 
compared  to  1,3,5-TNB  (Reference  53).  When  a  second  or  third  -Cl 
group  is  added  to  the  TNB  backbone,  there  is  no  further  effect  on 
the  rate  of  decomposition  (Reference  53) .  Substituents  with 
built-in  reactivity  (Reference  230)  can  override  all  of  the 
conventional  notions  about  thermal  stability  of  the  nitroaromatic 
ring.  Furthermore,  except  for  the  "ortho  effects"  discussed  in 
Section  IV,  the  MS  cracking  patterns  of  substituted 
trinitrobenzene  compounds  follow  few  systematic  patterns 
(References  126,170). 


SECTION  VI 


LATER  STAGE  REACTIONS 

All  of  the  forementioned  thermolysis  reactions  of 
polynitrobenzene  explosives  have  occurred  among  the  substituents 
while  leaving  the  aromatic  backbone  intact.  These  substituent 
homolysis  and  cyclization  reactions  are  obviously  early  reactions 
of  nitroaromatic  explosives,  and  can  be  thought  of  as  initiation 
steps  in  which  little  energy  is  generated.  The  next  stage  of 
thermal  decomposition  involves  propagation  reactions.  These 
reactions  involve  exothermic  oxidation  and  reduction  that  lead  to 
lower  molecular  weight  and  thermodynamically  stable  gaseous 
products.  While  much  is  known  about  the  related  kinetics  and 
thermochemistry  of  hydrocarbon/air  flames  (References  231-233)  and 
nitrogen  chemistry  in  combustion  (Reference  234),  the  processes  by 
which  large,  bulk  phase  energetic  molecules  thermolyze  to  form  the 
reactants  for  the  high  temperature  flame  is  a  challenging  area  of 
research  (References  41,107,217).  Specifically,  by  what 
intervening  steps  does  2,4-DNAn  or  the  dinitrotoluene  radical 
proceed  to  CO,  NO,  C02,  HCN,  N2  and  H20? 

At  least  in  the  gas  phase  according  to  MS/MS/CID  (Reference 
134)  ring  dissociation  does  not  occur  until  most  of  the  attached 
substituents  have  been  removed.  On  the  other  hand,  dissociation 
of  N02  from  the  ring  leaves  the  ring  with  excess  internal  energy 
to  drive  further  reactions  (Reference  111) .  In  the  simplest 
situation,  the  El  mass  spectrum  gives  a  scheme  for  further 
unimolecular  reactions  of  the  aromatic  ring  in  the  collisionless 
gas  phase.  2-Nitrotoluene  is  believed  to  react  further  according 
to  equation  (43)  based  on  2H  and  13C  labelling  (Reference  36) . 


Consistent  with  this  scheme  is  the  fact  that  anthranil  (equation 
12)  decomposes  by  rupture  of  the  NO  bond  with  a  rate  constant 
(Reference  34)  of  k  *=  3.7  x  1015  exp(-51084/T)  sec"1.  The  large  A 
factor  is  in  line  with  high  entropy  of  the  transition  state.  The 
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proposed  pathway  for  anthranil  degradation  is  equation  (44) 
(Reference  34) . 


The  relevance  of  the  findings  described  above  to  the  bulk 
state  of  an  explosive  is  dubious  because  of  the  important  role 
played  by  bimolecular  processes.  In  fact,  attack  on  the  ring  is 
quite  plausible  because  the  aromaticity  of  the  benzene  ring  is 


reduced  by  ring  closure  to  anthranil  (equation  12) ,  or 
furazan/furoxan  species  (equation  14).  This  makes  the  ring  more 
susceptible  to  reactions  that  ultimately  lead  to  C-C  bond 
cleavage.  For  example,  trinitroaromatics  are  notoriously 
susceptable  to  nucleophilic  addition  (References  235-237)  which 
also  reduces  the  aromaticity.  While  N02(g)  is  not  liberated  to 
any  great  extent  from  decomposition  of  trinitroaromatic  explosives 
at  temperatures  below  about  700 *C,  N02  or  NO  from  thermal 
decomposition  in  the  bulk  phase  may  be  attacking  the  aromatic  ring 
producing  CO  and  C02.  Maksimov  (Reference  63)  examined  the 
decomposition  of  TNB  vapor  and  noticed  first  a  rise  and  then  a 
drop  in  the  NO  concentration,  which  was  attributed  to  the  reaction 
of  a  nitrosoaromatic  compound,  RNO,  with  NO  according  to  equation 
(45) .  The  NO  consumption  occurs  as  C02  is  formed,  perhaps  owing 
to  the  reaction  of  the  nitrate  radical  with  *R  (Reference  63). 

RNO  +  2N0  - *  R-  +  N2  +  N03‘  (45) 

The  energy  released  by  the  propagations  reactions  of 
polynitroaromatic  explosives  is  sometime  reflected  in  the  relative 
concentrations  of  the  gaseous  products  that  are  liberated  after 
the  exotherm.  For  example,  Table  13  gives  the  IR  active  products 
for  TNB,  MATB,  DATB  and  TATB  (Reference  198) .  Endothermic 
products  NO,  HCN  and  N20  form  when  -NH2  groups  are  present  on  the 
ring.  HNCO  also  increases  as  the  NH2  content  increases.  Even 
though  the  C/O  ratio  in  this  series  of  parent  compounds  remains 
fixed,  the  COj/CO  ratio  increases  as  the  NH2  content  increases. 

The  energy  distribution  among  the  gaseous  products  reflects  the 
fact  that  the  heat  of  explosion  of  TNB  is  higher  than  those  of  the 
aminonitrobenzene  explosives  (Table  1) .  The  high  heat  of 
explosion  of  TNB  results  from  the  dominance  of  products  with 
negative  Ah;  (CO,  C02,  HNCO)  and  the  absence  of  products  with 
positive  AH;  (NO,  HCN,  N20) .  For  MATB,  DATB  and  TATB  the  amount 
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of  energy  from  the  exothermic  products  (CO,  C02)  is  lower  than  for 
TNB  and  is  further  aggravated  by  the  formation  of  endothermic 
products  (NO,  HCN,  N20) .  In  sum,  the  energy  output  of  TNB  results 
from  a  different  set  of  propagation  reactions  than  occurs  for 
MATB,  DATB  and  TATB. 
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Table  13.  Quantified  6aa  Products  from  Hasting  at  2000* C/sac  to 
520 *C  under  10-40  ata  Ar 

Percentage 


Coraomid 

C02 

CO 

NO 

HCN 

N20 

HNCO 

THB 

6 

94 

- 

- 

- 

- 

MATB 

17 

38 

28 

13 

- 

3 

DATB 

20 

31 

21 

18 

2 

8 

TATB 

27 

20 

26 

13 

3 

11 
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SECTION  VII 


CONCLUSIONS 

Consideration  of  the  literature  and  the  addition  of  new  data 
from  this  program  has  clarified  a  number  of  issues  about 
nitroaromatic  explosives.  The  main  points  are  highlighted  here. 

A  relatively  fruitless  approach  to  understanding  the  shock 
and  impact  sensitivity  of  nitroaromatic  explosives  is  the 
comparison  of  these  properties  to  the  parent  ground  state 
molecular  and  electronic  structure.  The  parameters  that  control 
shock  and  impact  sensitivity  in  this  type  of  compound  only 
partially  depend  on  the  molecular  and  electronic  structure  of  the 
compound.  Physical  properties,  such  as  the  void  volume  and 
density,  are  very  important.  The  essential  point  of  this  report 
is  that  the  rates  of  chemical  reactions  are  equally  important. 

The  analysis  presented  in  this  report  makes  it  completely 
understandable  why  the  time-to-explosion  does  not  correlate  with 
the  molecular  <;nd  electronic  characteristics,  such  as  are 
reflected  in  the  Hamir  >tt  or  Taft  linear  free  energy  relationships. 
The  time-to-explosion  at  times  longer  than  about  0.1  sec  is 
controlled  by  catalysis  of  the  decomposition  products  of  the 
parent  molecule  and  not  the  parent  molecule  per  se. 

In  general,  the  thermal  decomposition  of  nitroaromatic 
explosives  begins  with  one  or  more  endothermic  initiation  steps. 
These  steps  include  C-N02  homolysis,  isomerization  of  the  nitro 
group  to  the  nitrite  linkage,  and  reactions  of  the  non-energetic 
substituent  on  the  ring.  The  relative  importance  of  these 
initiation  steps  depends  on  the  temperature  and  pressure  reached 
during  the  stimulus  step.  Shock  initiation  predominantly  produces 
C-N02  homolysis  in  TNT  and  the  aminonitrobenzene  series.  Impact 
initiation  results  in  competitive  contributions  of  C-N02  homolysis 
and  reactions  involving  the  nonenergetic  substituent,  especially 
in  the  case  of  TNT.  In  sum,  decomposition  at  lower  temperature 
favors  reactions  of  the  non-energetic  substituent,  whereas  higher 
temperature  favors  direct  bond  homolysis. 

Following  the  initiation  reaction  steps  mentioned  above,  a 
different  gas/solid  product  mixture  will  exist  depending  on  the 
temperature.  The  propagation  steps  that  follow  will,  therefore, 
differ  and,  in  turn,  cause  the  explosive  characteristics  to  differ 
from  compound  to  compound.  For  example,  as  the  number  of  -NH2 
groups  increases  in  the  aminonitrobenzene  series,  the  quantity  of 
products  with  positive  heats  of  formation  increases  relative  to 
the  products  with  negative  heats  of  formation.  This  difference 
causes  a  decrease  in  the  heat  of  explosion  as  -NH2  groups  are 
added  to  the  ring. 
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It  would  appear  that  the  rational  design  of  a  nitroaromatic 
explosive  can  not  be  restricted  to  considerations  of  the  ground 
state  molecular  and  electron  structure  of  the  target  molecule. 
Instead,  the  initiation  and  propagation  reactions  need  to  be  the 
focal  issue. 
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The  Appendix 


The  appendix  contains  time-to-explosion  data  for  a  series  of 
nitroaromatic  compounds.  The  data  were  acquired  by  the  use  of  T- 
jump/FTIR  spectroscopy  (Reference  22).  About  200  ng  of  thinly 
spread  sample  was  heated  at  2000 *C/ sec  under  the  pressure  given  in 
the  caption.  The  time-to-exotherm  was  recorded  along  with  the 
gaseous  products,  which  were  quantified  by  the  well  established 
procedure  using  the  absolute  intensities  (Reference  41) .  Some  of 
these  data  were  used  in  the  text  of  this  report,  but  most  are 
simply  archived  here  because  they  are  not  discussed  but  are  valid 
time-to-explosion  data  in  the  time  range  given.  For  reasons 
discussed  in  the  text,  the  connection  between  the  time-to- 
explosion  and  the  parent  molecular  structure  is  vague.  The  time- 
to-explosion  is  predominantly  controlled  by  the  decomposition 
intermediates  of  these  compounds  which  in  turn  catalyze  the 
decomposition  further. 
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Figure  6.  T-Jump/FTIR  Data  for  1 , 3 , S-Triamino-2 , 4 , 6-trinitrobenzene  at  422 *C  and 
3.3  Atm  Ar.  Trace  CO,  appears  at  2.6  Seconds  and  Trace  H,0  was  Present 


Figure  7.  T-Jump/FTIR  Data  for  3 , 5-Diamino-2 , 4 , 6-trinitrotoluene  at  471 *C  and  3.3  Ata  Ar 
Trace  C02  appears  at  0.6  Seconds  and  Trace  H20  was  Present. 
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Figure  9.  T-Jump/FTIR  Data  for  2 , 4 , 6-Trinitroresorcinol  at  521 *C  and  3.3  Atm 
H,o  was  Present. 
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T-Jump/FTIR  Data  for  1, 3, 5-Triamino-2, 4, 6-trinitrobenzene  at  390*C  and 
6.67  Atm  Ar.  Trace  H,0  was  Present. 
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Figure  12.  T-Jump/FTIR  Data  for  2,4,6-Trinltroaniline  at  454 *C  and  10  Atm  Ar 
Trace  H,0  was  Present. 
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Figure  14.  T-Jump/FTIR  Data  for  3 -Amino-2, 4, 6-trinitrotoluene  at>513‘C  and  10  Atm  Ar 
Trace  H,0  was  Present. 
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T-Jimp/PTIR  Data  for  1, 3, 5-Triaraino-2 , 4, 6-trinitrobenzene  at  452 
and  10  Atm  Ar.  Trace  H20  was  Present. 
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Figure  18.  T-Jump/FTIR  Data  for  2,3,4,6-Tetranitroaniline  at  308  *C  and  10  Atn  Ar 
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Figure  21.  T-Jump/FTIR  Data  for  Picrylhydrazine  at  323  *C  at  10  Atm  Ar.  KNCO  appears  at 
0.80  Seconds  and  Trace  H,0  was  Present. 
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Figure  22.  T-Jump/FTIR  Data  for  Picric  Acid  at  412 *C  and  20  Atm  Ar 


Figura  23.  T-Jump/FTIR  Data  for  Picric  Acid  Potassium  Salt  at  587 *C  and  20  Atm  Ar 
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Figure  25.  T-Jump/FTIR  Data  for  1,3 -Diamino-2, 4, 6-trinitroben*ene  at  375 *C  and 
20  Atm  Ar.  Trace  Hz0  was  Present. 


Figure  26.  T-Jump/FTIR  Data  for  1, 3, 5-Triamino-2 , 4 ,6-trinitrobenzene  at  402 *C  and 
20  Atm  Ar.  Trace  H,0  was  Present. 
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Figure  27.  T-Jump/FTIR  Data  for  2,4,6-Trinitrotoluene  at  405*C  and 
20  Atm  Ar.  Trace  CO,  appears  at  1.1  Seconds. 
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Figure  29.  T-Jump/FTIR  Data  for  2-Picrylpropane  at  484 *C  and  20  Ate 
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Figure  30.  T-Jump/FTIR  Data  for  2 ,4,6-Trinitro-£-butylbenzene  at  405 *C  and  20  Ate  Ar 
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Figure  31.  T-Jurap/FTIR  Data  for  Picrylhydrazine  at  300 *c  and  20  Atm  Ar.  Trace  H20  was 


Figure  32.  T-Jurap/FTIR  Data  for  2 , 4 , 6-Trinitroanisole  at  326*C  and  20  Ata  Ar 
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Figure  33.  T-Jump/FTIR  Data  for  2 , 4 , 6-Trinitrobenzamide  at  510 *C  and  20  Atm  Ar 
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Figure  34.  T-Jump/FTIR  Data  for  2 , 4 , 6-Trinitroacetophenone  at  506 *C  and  20  Ata  Ar 
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Figure  35.  T-Ju*p/FTIR  Data  for  Methyl  2,4, 6-Trinitrobenzoate  at  481*C  and  40  Ata  Ar 
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,gure  36.  T-Jump/FTIR  Data  for  1,3,5-Trinitrobenzene  at  470*C  and  40  Atm  Ar 
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Figure  37.  Tine-to-Exotherro  Data  and  Apparent  Activation  Energy  of  2,4, 6-Trinitroaniline 
at  10  Atm  Ar. 
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Figure  38.  Time -to-Exo therm  Data  and  Apparent  Activation  Energy  of  2,3,4, 6-Tetranitro- 
aniline  at  10  Atm  Ar. 


(XJLX)ttl 


119 


Figure  39.  Time-to-Exotherm  Data  and  Apparent  Activation  Energy  of  2,4, 6-TrinitrobenEamide 
at  10  Atm  Ar. 
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Figure  40.  Tlme-to-E*>ther»  Date  end  Apparent  Activation  Energy  of  Picrylhydr.aine  at 
10  Atm  Ar. 
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Figure  41.  Time-to-Exotherra  Data  and  Apparent  Activation  Energy  of  Picric  Acid  at 
10  Atm  Ar. 
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Figure  42.  Time-to-Exotherm  Data  and  Apparent  Activation  Energy  of  2 , 4 , 6-Trinitroaniline 
at  20  Atm  Ar. 
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Figure  43.  Time-to-Exotherra  Data  and  Apparent  Activation  Energy  of  2, 4, 6 -Trinitrotoluene 
at  20  Atm  Ar. 
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Time-to-Exotherm  Data  and  Apparent  Activation  Energy  of 
2 , 4, 6-Trinitroethylbonzene  at  20  Atm  Ar. 
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Figure  45.  Tirae-to-Exotherm  Data  and  Apparent  Activation  Energy  of  2-Picrylpropane 
at  20  Atm  Ar. 
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Figure  47.  Time-to-Exotherm  Data  and  Apparent  Activation  Energy  of  Picrylhydrazine 
at  20  Atm  Ar. 
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LOW  TEMP.  ACTIVATION  ENERGY 
=  14.7  KCAL/MOLE 


:o-Exotherm  Data  and  Apparent  Activation  Energy  of 
■Trinitroacetophenone  at  20  Atm  Ar. 
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Figure  51.  Tiae-to-Exothera  Data  and  Apparent  Activation  Energy  of 
Hethyl  2 , 4 , 6-Trinitro-£-butylbenzene  at  20  Ata  Ar. 
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,  5-Trinitrobenzene  at  40  Atm  Ar. 
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